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Abstract. These lectures stress the theoretical elements that underlie a wide 
range of phenomenological studies of high-energy QCD, which include both soft 
and hard processes. After a brief introduction to the basics of QCD, various 
aspects of QCD-based phenomenology are covered: colour transparency, hadronization 
of colour charges, Regge phenomenology, parton model, Bjorken scaling and its 
violation, DGLAP evolution equation, BFKL formalism, GLR-MQ evolution equation 
and saturation. In the last part of the lecture, we employ the light-cone dipole 
formalism to describe deep inelastic lepton scattering, Drell-Yan processes, direct 
photon production, diffraction, quark and gluon shadowing in nuclei, the Cronin effect 
and nuclear broadening. 
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1. The theory of strong interaction 



Strong interactions are described by a quantum field theory known as quantum 
chromodynamics (QCD). In many ways QCD is a unique theory. Quantum 
electrodynamics (QED), and its expansion to the electroweak Standard Model of particle 
physics, is also a quantum field theory. QED is a renormalizable theory but it loses all 
its credibility as we approach the energy scale, the so-called Landau pole position, where 
the strength of the coupling constant (the strength of the interaction) becomes infinit^. 
On the other hand, if the cutoff goes to infinity, QED becomes trivial. QED is not the 
only theory with a Landau pole problem; every theory which is not asymptotically 
free suffers from this problem. QCD is the only known theory which is free from such 
problems. QCD needs only a few parameters to be defined completely: one universal 
coupling strength and one mass for each kind of quark. 

Despite more than half a century of attempts, our knowledge about many aspects of 
QCD is still rudimentary. This is mainly due to the fact that QCD evolves from a few- 
body theory of free quarks and gluons at short distances to an extremely complicated 
infinite-body theory of objects like hadrons and nuclei, giving rise to a variety of complex 
physical systems and their interactions. The aim of this manuscript is to bring together 
various aspects of high-energy nuclear physics as tools for studying QCD itself. 



1.1. QCD Lagrangian and it symmetries 

The Lagrangian of QCD is given by 

C = q{iYd, - m')q - \{F;,f + gq^A^q, (1) 

where q is the quark field which is defined in the fundamental representation of the 
colour and fiavor group, and the conjugate Dirac field is defined as g = g^7°. The gluon 
field matrix A^ = y4^A"/2 is defined in the fundamental SU{Nc = 3) representation 
where Nc denotes the number of colour. A" being the generators of the gauge group 



I The dependence of coupling constants on the energy scale is one of the basic ideas behind the 
renormalization group which will be discussed in Section 1.3. 
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Figure 1. Gluons carry colour charge and interact witli each other via these vertices. 

which satisfies [A"/2, = if"-^^\^/2 where f"-^^ are the structure constants of S't/(3). 
We define g as the strong couphng constant. The field strength F^^ is given by 

F;, = d,Al - d^Al + gr'^AlAl. (2) 

The non-Abehan nature of QCD is manifested by the quadratic term in the gauge 
field strength, which gives rise to gluon-gluon interactions shown in Fig. [H The crucial 
difference between QCD and QED is the presence of this quadratic term which makes the 
QCD field equations non-linear. These nonlinearities give rise to a non-trivial dynamics 
and various rich structures which are unique properties of the strong interaction. The 
colour and fiavour indices of the quark field are suppressed. m° is the current quark 
mass which is not directly observable if QCD confines quarks. The current quark mass is 
colour independent and can be brought diagonal in fiavour space. There are six fiavours 
of quarks, each of which has a different mass. The three light quarks are called up (u), 
down (d) and strange (s), while the three heavy quarks are called charm (c), bottom 
(b) and top (t). The following values for the light current quark masses are found in 
the Particle Data tables p^J, 

ml = 2 to 8 MeV, = 5 to 15 MeV, m° = 100 to 300 MeV. (3) 

Notice that the quark masses are renormalization-scheme dependent. The above 
values are obtained in a subtraction scheme at a renormalization scale O(lGeV). In 
addition to flavour, quarks carry another quantum number known as colour. Each 
quark comes in three colours which, based on a convention, are called red, green and 
blue. 

The Lagrangian Eq. ([1]) has a large classical symmetry: we have the local gauge 
symmetry SU (iVc) by construction, 

q ^ Us. q ^ qUl U,{x) = exp(zr 

A^ ^ U.A^Ul - -U.id^Ul 

In QED, there is only one electric charge, and the gauge transformation involves a single 
phase factor U = exp{ia{x)). The QCD Lagrangian Eq. ([T]) has also a global flavour 
symmetry which does not affect the gluon flelds, 

q^Uvq, q^qUl, [/y = exp(^^^(^)^). (4) 
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where (^)f denotes the generators of the flavour group U{Nf) and Nj denotes the 
number of flavors. The above symmetry is referred to as vector flavor symmetry 
Uv{Nf). When the generator is the unit matrix, we have Uv{^) symmetry associated 
with conservation of baryon number. There is another global symmetry which is exact 
at m° = 0, namely chiral symmetry. This symmetry is very similar to vector flavor 
symmetry, apart from an extra factor of 75 in the generator of the transformation. 

q UAq, q ^ qUA, Ua = exp ^^Ts^'^ly • (5) 

Notice that due to the factor 75 the quark fleld and its conjugate partner are transformed 
by the same matrix in contrast to vector transformation Eq. (jl]). This transformation 
Eq. is called the axial-vector transformation and can be combined with the vector 
transformation to deflne a bigger symmetry at chiral m° = which is then called chiral 
symmetry Uv{Nf) x UA{Nf). One may alternatively deflne right- and left-handed quark 
flelds by following transformation 

1-75 1 + 75 

The right- and left-handed massless fermions are eigenvalues of the helicity or chirality 
(with eigenvalue ±1) and are not mixed together. The chiral symmetry can be 
equivalently written as Ul^Nj) x Ur^Nj). 

Not all the above-mentioned symmetries survive quantization. Particles with 
opposite helicity are related by a parity transformation, therefore in a chirally symmetric 
world, the hadrons should come in parity doublets. However, in real life we do not 
observe such degeneracy. Therefore one can conclude that chiral symmetry is not 
realized in the ground state and chiral symmetry is spontaneously broken. A theory 
where the vacuum has less symmetry than the Lagrangian is called a theory with 
spontaneous symmetry breaking. The Goldstone theorem [6] tell us that the spontaneous 
breaking of a continuous global symmetry implies the existence of associated massless 
spinless particles. This indeed was conflrmed due to the existence of the light 
pseudoscalar mesons in nature (pions, kaons and etas) which may be assigned as 
pseudo-Goldstone bosons [2]. Moreover, the existence of a quark condensate {qq) 
implies that the SU^Nf)^ x SU{Nf)fi symmetry is spontaneously broken down to 
SU{Nf)v- Therefore one may conceive QCD quark condensate as an order parameter 
for chiral symmetry breaking. The concept of spontaneous broken chiral symmetry is 
the cornerstone in the understanding of the low-energy hadronic spectrum. 

The U{1)a symmetry implies that all hadrons should come with opposite parity 
partners. However, this is not the case, therefore this symmetry must be broken 
somehow. If the spontaneous symmetry breaking mechanism works here, then one 
should observe a Goldstone boson associated with ^7(1)^, namely an / = pseudoscalar 
meson having roughly the same mass as the pion. Surprisingly there is no such Goldstone 
boson. This problem is sometime called U{1)a puzzle. It turned out that the U{1)a 
symmetry is explicitly broken by quantum effects. This effect is known as the axial 
anomaly [3j. The axial charge corresponding to the axial current = q'jfj.'y^q is not 
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Figure 2. The diagram corresponding to the C/4(l)-anomaly. 

conserved because of the contribution of the triangle graph in Fig.[2l The four- divergence 
of the axial current is given by 

Nt 

d^Jl = Y: 2im,q^'q + ^trG^^G^,, (7) 

where G^j, = e^j^^kxC^^ /2 is the dual field strength tensor. The last term (gluonic part) 
is a full divergence, and one may expect that this term not to have any physical effect 
if the QCD vacuum were trivial. It was shown by 't Hoof that due to instanton effects, 
the U{1)a symmetry is not manifested in nature |3j. 

Finally, at m° = 0, the QCD Lagrangian is invariant under a scale transformation 
which is called dilatational symmetry: 

q{x) ^ t''\{t'^x), Al{x)-^tA';^{t~^x), x,-^t-^x,. (8) 

This symmetry is again broken at the quantum level due to the trace anomaly [4]. 

1.2. QCD versus QED 

Let us remember the main differences between QCD and QED. QCD is an extended 
version of QED which now, instead of one charge, has three different kinds of charge 
called colour. Similar to the photon in QED, here massless spin-one particles, the gluons, 
respond to the presence of colour charge. The colour charged quarks emit and absorb 
gluons in the same way as electrically charged leptons do. However, radiation of a photon 
does not change the charge of the electron, while a gluon can change the quark colour. 
The response of gluons to colour charge, as measured by the QCD coupling constant, 
is much more drastic than the response of photons to electric charge. Gluons, unlike 
photons, interact directly with each other, although the colour charges, like electric 
charge in QED, are conserved in all physical processes. Therefore gluons must be able 
to carry unbalanced colour charges in contrast to their counterpart the photon in QED. 

In the following sections we shall also recapitulate the most important features 
of QCD which are not accessible perturbatively. These non-perturbative features are 
unique for QCD and should be traced back to the main differences between QCD and 
QED. 
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Figure 3. Right: The running couphng constant as a function of momentum transfer 
determined from different processes. Left: Summary of ag [9]. 



1.3. Asymptotic Freedom 

Having introduced the gauge fixing term and an associated ghost term by means of 
the Faddeev- Popov procedure [5], E], one can carry out perturbation theory in terms 
of coupling. Similar to QED, a dimensionless physical quantities TZ can be expressed 
by a perturbation series in powers of the coupling parameter as {as is the notation 
for (yf^/47r). Owing to the renormalization process, a renormalization scale /i enters the 
algebra [7j in order to remove the ultraviolet divergence. Therefore, one can write the 
dimensionless quantities TZ in terms of other available dimensionless parameters Q^/iJ^'^ 
and the renormalized coupling ^^(/i^). However, the physical quantity TZ cannot depend 
on the arbitrary /x. This means that TZ should be renormalization scale invariant 



dU 



d oda^ dd 



This equation exphcitly shows that any dependence of 7^ on /i must be cancelled by an 
appropriate /i-dependence of a^. It is also natural to identify the renormalization scale 
with the physical energy scale of the process, i.e. fi"^ = Q^. The running coupling is 
described by the renormalization group equation [7], 

Q'^=p{»sm) ■ (10) 

Whenever the coupling is small, the (3 function can be computed perturbatively, 

f3{as) = -^aim - ^,aim + (11) 

with 

_ 33 - 2Nf _ 153 - 19iV, 

- 127r ' - 247r2 " ^^^^ 
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Therefore one can readily calculate the effective running coupling at one-loop level 
ignoring the /3i term, 

as{Q') = — (13) 

where A is a scale parameter of QCD and depends on the subtraction scheme and 
the number of active flavours. The present world average for at the mass is 
as{Mz) = 0.118 ± 0.002 which leads to 

= (208+i)MeV, (14) 



where the symbol MS stands for minimal subtraction scheme [7] and the superscript 
indicates the number of active flavours. This value is taken from an analysis of various 
high energy processes [8[ [9] , see also Fig. [31 The most striking feature of the running 
coupling is that it decreases logarithmically with for Nf < 17 when /3o > 0. This 
originates from the self-interaction of gluons which leads to anti-screening, in contrast to 
QED where the sign of /3o is negative. Therefore perturbation theory works very well for 
large Q^. This phenomenon is called asymptotic freedom [10]. However, if is near 
^M5' perturbation theory does not work anymore and non-perturbative phenomena 
enter the stage. One of the biggest challenges of QCD is to connect these two domains. 
Admittedly, there is yet no unambiguous method to connect small and large distances 
in QCD. 



1.4- Chiral symmetry breaking 

In the first section, we introduced the symmetries of the QCD Lagrangian. In the limit 
of massless quarks, QCD possesses chiral symmetry UiiNf) x Uii{Nf) which means that 
left- and right-handed quarks are not mixed, 

qL^VLQi; qR^VnqR- VL,VReU{Nf). (15) 

As we already discussed, owing to the presence of the quark condensate (gg), chiral 
symmetry is spontaneously broken and left- and right-handed quarks and antiquarks 
can transform into each other: 

{qq) = {qiqu) + {quqi)- (16) 

Dynamical chiral symmetry breaking is one of the important non-perturbative features 
of QCD which is responsible for the generation of quark massei§|. In order to show that 
this phenomenon is purely non-perturbative, we employ the QCD gap equation [T3|, 

S{p)-^ = {il.p + m') + I ^yD,^{p-q)^j,S{q)T:{p,q), (17) 

§ There is another very different way to generate mass from vacuum, the so-called Casimir effect 
[llj , which originates from the response of the vacuum in the presence of non-perturbative boundary 
conditions. The existence of boundary conditions in quantum field theory is not always free of problems 
(see, for example, Ref. [12]). 
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where imP and g are the current-quark bare mass and the couphng constant, respectively. 
D^y{p—q) is the dressed-gluon propagator and r"(p, q) is the dressed-quark-gluon vertex. 
The general solution of the gap equation is a dressed-quark propagator of the form 

o,. 1 _ Z{p^) 

i-f.pA{p^) + B{p^) i^.p + M{p^)' ^ ' 

The functions ^(p^) and -B(p^) contain the effects of vector and scalar quark-dressing 
induced by the quark interaction with the gluon field. The function M(p^) denotes the 
quark mass. One may now use the gap equation to work out the fermion self-energy 
perturbatively [Tl|. One obtains, 

= ^0 _ ^ ln(pVm2) + ...^ . (19) 

It is observed that at all orders of the loop expansion, terms are proportional to the 
current-quark mass and consequently vanish as m° 0. The quark mass is defined 
as a pole of the dressed-quark propagator; therefore no mass is generated at a current- 
quark mass equal to zero, i.e., the dynamical chiral symmetry breaking is impossible in 
perturbation theory and there is no mixing between left- and right-handed quarks at the 
perturbative level. Notice that, apart from the trivial solution -B(p^) = at m = 0, a 
non-trivial solution -B(p^) 7^ can indeed be found at the chiral point, albeit accessible 
non-perturbatively. The renormalization effect is not included in Eq. ( |T7|) . but it does 
not change the above argument [H]. The quark condensat^ in QCD is given by the 
trace of the full quark propagator Eq. 



(gg) = —i lim Tr S{x, y). (20) 
y 

Notice that since qq is a gauge invariant object, one may take any gauge to obtain 
the dressed quark propagator which has a general form as equation ( fTSj) . It is obvious 
that when -B(p^) = 0, the quark condensate does not take place, simply because of the 
identity Tr7^ = 0. It has been shown in many non-perturbative approaches that the 
emergence of a dynamical quark mass leads to the non- vanishing of quark condensate 
and vice versa, see, for example, Refs. [T7t 118]. 



1.5. Confinement 

Another important non-perturbative feature of QCD is colour confinement [10] . Loosely 
speaking, confinement is defined as the absence of any free coloured objects in nature. 
But it is possible that there exists a composite coloured particle which can form 
colourless bound states with another coloured particle like quarks. Colour confinement 
is still not properly understood, and a clear and indisputable mechanism responsible 
for this effect remains yet to be discovered. The basic property of confinement can be 

II Note that, at finite density and temperature, the formation of a quark cooper pair condensate 
(gq) 7^ is also possible, leading to colour symmetry breaking, the so-called colour superconductivity 
phenomenon(BCS) [IS] and diquark Bose-Einstein condensation(BEC) ^16) . 



Applied High Energy QCD 



10 



explored by looking at heavy qq propagation at a large distance i? in a time interval T. 
The behaviour of such a system can be described by the Wilson loop, 

W{R,T) =Tr[Pexp (^i J^A^^T'^dx^'^], (21) 

where T*^ denotes the generator of SU{3). One can show that at large interval of time 
T, 

W{R,T ^oo) =exp {-TV (R)), (22) 

where V{R) is the static potential between the heavy quarks. At large distances this 
potential grows linearly: 

V{R ^ oo) = aR. (23) 

Therefore the Wilson loop at large R and T behaves as W{R oo,T oo) = 
exp {—aTR), which is the so-called area low and indicates confinement. 

Confinement originates non-perturbatively, since it is associated with a linear 
potential with a string tension 

a(xA^exp{- [ ) , (24) 



m, 

which is obviously non-perturbative in the coupling. Note that the string picture of 
quark confinement is not free of flaws, since string breaking will occur once the potential 
energy approaches the quark pair creation threshold. 

It is well-known that for confinement it is sufficient that no coloured Schwinger 
function possesses a spectral representation. This is equivalent to say that all coloured 
Schwinger functions violate reflection positivity [13] . Another way of realization of QCD 
confinement is due to Gribov theory in which colour confinement is determined by the 
existence of very light (almost massless) quarks [19]. There are in fact many different 
ways that the confinement can be realized, such as monopole condensation, infrared 
enhancement of the ghost propagator, etc. For a review of this subject see Ref. [20] . 

One may wonder if there is a non-trivial solution for the gap equation B{p^) ^ 
which gives rise to a pole of the quark propagator, this might contradict QCD 
confinement since the quark is coloured. Indeed this is one of the subtle point in 
every QCD model and cannot be easily resolved. In principle, there will be a long- 
range force between massive quarks to confine them and also a short range spin-spin 
interaction between massive dressed quarks. The former will modify the low momentum 
part of the propagator to remove the quark from being on-shell. Actually, this describes 
a phenomenologically motivated picture of a constituent quark model based on the 
dynamical symmetry breaking. Having said that, it is very hard to incorporate the 
dynamical symmetry breaking and the confinement into a QCD model. In fact, many 
models constructed to describe the low-energy properties of hadrons [18], El] are assumed 
to be only dominated by the quark flavor dynamics and dynamical symmetry breaking 
and are indeed reliable only at intermediate scales, between confinement scale few 
hundred MeV up to a scale about 1 GeV. 
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Figure 4. e+e annihilation to qq or pair. 



2. Evidences for coloured quarks 

Historically, the idea of colour degree of freedom emerged as a viable solution to the 
problem of how to construct the wave function for the doubly charged baryons |22j . 
The wave function of in space, spin and flavour is symmetric and violates the Pauli 
exclusion principle since A"*""*" is fermion with spin 3/2. This problem was resolved by 
introducing a new degree of freedom, the colour degree of freedom, and requiring that 
the A^^ wave function to be antisymmetric in the colour degree of freedom. 

Although coloured states are not detected in experiments, and only colour singlet 
states exist in nature, there is much experimental evidence in favour of a colour degree 
of freedom. One of the direct experimental test for a colour degree of freedom comes 
from e'^e~ annihilation into hadrons. In the e^e~ annihilation process, first a pair of 
quarks e^e^ qq is produced which then fragment into hadrons. The cross-section for 
producing a free qq pair is the same as for producing a fi^fi~ pair except for the quark 
charge and colour number which should be replaced with the muon charge, see Fig. HI 
Therefore in order to extract information about the QCD content of e~^e~ annihilation, 
in particular, the colour degree of freedom, it is convenient to express the total cross- 
section of e~^e~ —>■ qq annihilation in units of the cross-section of n production, 

e~^e~ — > Hadrons 

R = I _ • (25) 

e+e — i> 

The cross-section to produce any number of hadrons is proportional to that to 
produce a fi~^fi~ pair. This is because a highly virtual photons decays to quarks in a time 
scale t ~ (where a/s is the center of mass energy), while a hadron with mass 

needs a formation time t ~ l/M^. Therefore, there is not enough time for confinement 
to affect the annihilation cross-section and one can assume that the produced qq pair 
fragments into hadrons with unit probability. 



(26) 
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Figure 5. R in the light-flavor, charm, and beauty threshold regions. Data errors 
are total below 2 GeV and statistical above 2 GeV. The full list of references to the 
original data and the details of the R ratio extraction from them can be found in [23j . 



Therefore one finds 



R = N, J2 



(27) 



-uA,. 



where the factor Nc is the number of colour and Cq denotes the quark charge. The 
summation in the above equation is over all flavours that are kinematically allowed. 
Depending on energy, various flavour degrees of freedom contribute, 

2 

-N^ (u, d, s). 



R 



9 

9 ' 



[u, d, s, c). 



[u, d, s, c,b). 



(28) 



up to about 3 GeV only u, d and s contribute, while at higher energies charm and b 
quarks start contributing as well. If one assumes that A^^c = 3, then Eq. ( l28l) predicts 
i? = 2, y and y, respectively. If we ignore nonperturbative effects close to threshold, 
such as the formation of bound states, we expect R to present a series of steps as a 
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Figure 6. Decay of pion to photons. 

function of ^/s. In Fig. [5] we show various experimental data which shows remarkable 
agreement with more detailed perturbative QCD calculation based on the assumption 
that = 3. 

Another strong evidence of colour degree of freedom is the measurement of the 
neutral pion decay into photons tt^ — > 77. The pion decay rate is computed from 
the triangle diagram shown in Fig. [6l Because of the quark loops, the decay rate is 
proportional to N^. The experimental value of the pion decay rate can only be described 
byiV, = 3 [21]. 

3. Colour transparency (CT) 

So far we have treated colour as just a new quantum number, a new degree of freedom. 
Is there any evidence that this colour is responsible for the strong interactions? 

If an interaction is controlled by colour, how can colourless hadrons interact? 
Apparently, only due to the spatial distribution of colour (carried by quarks and gluons) 
inside the hadrons, i.e., due to the existence of hadronic colour-dipole momentum. 

This observation immediately leads to experimentally observable consequences. 
Since colourless dipoles of vanishing size cannot interact, the interaction cross-section 
of such a dipole (say, quark-antiquark) with other hadrons should vanish when the 
transverse dipole separation goes to zero [25] . 

airr) oc . (29) 

This remarkable relation deserves commenting upon: (i) only transverse dipole 
separation matters, since at high energies longitudinal momentum transfer in exclusive 
reactions vanishes. For example, if the beam particle of mass mi and energy E (in the 
target rest frame) is excited to mass m2, while the target remains intact, the longitudinal 
momentum transfer reads, = {m\ — m\)/2E] (ii) the quadratic rr-dependence is 
dictated by dimension counting, no other dimension parameters can be used here (the 
QCD scale Aqcd may enter only via the coupling as); (iii) an additional logarithmic 
dependence on may and does exist [25]; (iv) such a small-rr behaviour is common 
for QED and QCD; however, in the former case the total cross-section is predominantly 
elastic, while in the latter case is inelastic. 





Figure 7. Top: decp-inclastic electron proton scattering, ep eX, at large Bjorken 
x: the virtual photon knocks a valence quark out of the proton, whose remnants form 
the final hadronic state X. Bottom: elastic ep — *■ e'p': the initial proton caught in a 
small size configuration survives a strong kick with increased probability. 



3.1. Quasielastic scattering off nuclei 

The experimentally measured total hadronic cross-sections is a result of the interplay of 
different dipole sizes whose probabilities are controlled by the hadronic wave functions. 
In some cases the probability of small size configurations in a hadron can be enhanced 
leading to a reduced interaction cross-section of such a hadron. An example is elastic 
electron-proton scattering, ep e'p', with high momentum transfer [26l [27], as 
illustrated in Fig, [3 

When the reeoil protou has a reduced sizfl it should interact more weaker than a 
regular proton with other targets. Such a possibility exists, if the elastic ep scattering 
is embedded into a nucleus, i.e. in quasielastic y4(e, e'p) A* reaction. The benchmark 
to compare with in this case is the expectations based on Glauber model calculations, 
where the recoil proton attenuates exponentially with the path length in the nucleus 
and with the normal proton-nucleon cross-section. 

The observable usually measured in such experiments is nuclear transparency, 
defined as, 

Z a[ep e'p') 

Basing on the above ideas of colour transparency one should expect a deviation rising 
with from the Glauber model predictions. Unfortunately no experiment performed 
so far has provided a clear evidence for such an effect. The results of the dedicated 
experiment NE18 at SLAG [2H] are depicted in Fig. [Sl^right). 

Apparently data show no preference either for Glauber, or GT based models |29j . 
Other measurements of A(e, e'p) A* reactions were not successful either, when searching 
for a GT signal. To fit the cross-section of this reaction by power A- dependence, 
(j(e, e'p) oc A", a rise of a with would be a signal of GT. However, the collection of 
data [30j depicted in Fig. [H](left) versus show no rise. Moreover, the value a = 0.75 
agrees with what one should expect from the Glauber model. 

^ Strictly speaking this is not a proton. This state can be projected either into a proton (as in the 
present case), or proton excitations. 
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Analogously, experiment in quasielastic proton-proton scattering, A(jp,2p)A*, 
performed at BNL |3T] did not provide any clear signal of CT. Although data deviates 
from the Glauber model predictions, at higher momentum transfers the agreement is 
restored. 
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Figure 8. Right: Nuclear transparency measured in quasielastic scattering on iron 
(upper panel) and gold (lower panel) in NE18 experiment at SLAC (28J. Dashed and 
solid curves present expectations based on the Glauber model and CT [29|. Left: 
Data from different experiments on quasielastic electron scattering, A{e,e'p)A* , for 
A-dependence of the cross-section fitted by A"' [30]. The curve is a = 0.75. 



Why did these experiments fail to observe a CT effect? It turns out that it is not 
enough to produce a small-sized configuration in a hard reaction. The produced hadron 
has to maintain this small size during propagation through the nucleus. It is clear that 
in a sufficiently long time interval the hadron will develop its wave function and restore 
the regular size. The time scale controlling this process is called formation time and for 
a recoil proton is given by, 

2E 

^/<-^^-0.4fmxE,(GeV). (31) 

Here Bp is the energy of the recoil proton , and rup* is the mass of the first proton 
excitation. In order to have // 3> Ra for heavy nuclei Ra ~ 5fm, the proton energy 
should be much higher than 10 GeV. The highest energy of recoil protons in the NE18 
experiment [28] was Ep ^ 4 GeV which is too low to keep the size of the produced 
hadron small within the nuclear range. 

This is the principal problem of quasielastic scattering where the photon energy 
u and virtuality are strongly correlated, 2mpV = Q^. Thus the recoil proton energy is 
Ep ^ Q'^/2mp. Therefore one must go to extremely high virtualities, > 20GeV^, 
just in order to increase Ep. However, the cross-section becomes vanishingly small. 
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3.2. Diffractive electroproduction of vector mesons 

Diffractive virtual photoproduction of vector mesons is free of this problem. The space- 
time development of this reaction is illustrated in Fig. [9] 



V 



Figure 9. Virtual diffractive photoproduction of vector mesons. A virtual photon 
fluctuates into a qq pair of transverse separation r"^ ~ which propagates through a 
nucleus, interacts diffractively, and being brought to the mass shell develops the wave 
function of the vector meson. 

At very high energies, the qq fluctuation lifetime, 

= Q^+Mj' ^^^^ 

(which is also called coherence time), becomes very long. So one can treat the qq dipole 
propagating through the nucleus as "frozen" by Lorentz time dilation at the initial size 
rf. ~ Thus, one can keep the scale finite, while the photon (and vector meson) 

energy can be increased with no restriction. This is the main advantage of this process 
for the search for CT effects compared to quasielastic reactions. The first measurements 
proposed in [32] and performed by the E665 collaboration [35| confirmed the theoretical 
expectations [33] of CT effects depicted in Fig. fTOT right). The high photon energy in this 
experiment led to Ic ^ Ra and allowed us to greatly simplify theoretical calculations. It 
turns out that at the opposite limiting case of 1^ <^ Ra and high Q^, the photon energy 
may be still high enough to keep the formation time scale Eq. (13T1) sufficiently long to 
observe CT effects. In this case a signal of CT would be a rising energy dependence 
of nuclear transparency. Corresponding measurements are under way at Jefferson Lab 

Similarly to diffraction, quasi-free hadron scattering off a nucleus can be performed 
at high energies, while the hadron size can be controlled by transverse momentum 
p^ —t [35]. In the case of Reggeon exchange, the pion formfactor formfactor 
suppresses large-sized configurations in the hadronic wave function at rather small t. 
Measurements were performed by the PROZA collaboration [36] with 40 GeV pions 
in quasi-free charge exchange scattering tt^A tt^A*. The results are depicted 
in Fig. [Tor left) in comparison with Glauber model expectations (dashed curve) and 
calculations including CT effects [35] . 

Notice that both models predict a peak at —t ^ 0.6 GeV^, because the cross- 
section of free scattering, 7r~p — > vr^ra, has a minimum at this momentum transfer, and 
the position of the minimum in quasi-free scattering is shifted by multiple interactions 
in the nucleus. 
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Figure 10. Right: Nuclear transparency as function of for carbon, calcium and 
lead. Data points from the E665 experiment at Fermilab [34j are compared with 
calculations [33]. Left: Nuclear transparency in quasi- free charge exchange of pions 
on carbon. Data are from [36]. Solid and dashed curves represent calculations [35] 
including or disregarding CT effects, respectively. 



4. Bags, strings... 

Gluonic condensate in vacuum puslies tlie energy density below tlie perturbative level, 
^vac < 0. If the colour field of the valence quarks suppresses vacuum fluctuations, then 
the energy density inside the hadron is higher than outside. Therefore the vacuum tries 
to squeeze the hadron. However, the chromo-electromagnetic energy [E"^ + H'^)/2 inside 
a smaller volume rises leading to an equilibrium, as illustrated in Fig. [TTT a). Thus 
hadrons look like bubbles in the QCD vacuum, this is the key idea of the the MIT bag 
model [37] . 

What happens if a quark is knocked out with a high momentum? On account of 
the same properties of the QCD vacuum the chromo-electric flux is squeezed into a tube 



evac<0 




t t t t t t t 1 1 t t t t T t ^ 



a 



Figure 11. a: Pictorial illustration for the MIT bag model; b: a stretched bag 
becomes a tube of a constant cross-section, which can be treated as a string. 
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of a constant cross-section, 

vrr^ = f , (33) 

as illustrated in Fig. [TlT b). Here g is the colour charge at the ends of the tube; k is 
the energy density stored in the tube per unit of length. This pattern of colour fields is 
quite different from that in QED, as illustrated in Fig. [121 



QED QCD 

V(r) = ^ V(r) = K L 




Figure 12. a: Electric field pattern in QED: the potential falls with charge separation 
as 1/ L\ b: Colour field pattern in QCD: the field is squeezed into a tube which breaks 
up by production of qq pairs tunnelling from vacuum. 



The potential between two electric charges falls with distance as 1/L, while in 
QCD it rises linearly. In fact, the rising with distance of a string potential explains the 
observed linearity of Regge trajectories (see below). 

Usually the transverse size is not important, so the tube may be treated as a one- 
dimensional string, and k is called string tension. It can be either calculated on the 
lattice, or related to the universal slope of Regge trajectories = 0.9 GeV~^ [38] . 

1 1^. (34) 



27ra^ fm 

This energy is sufficient for the creation of a couple of constituent quarks via tunnelling 
from the vacuum. One can hardly stretch a string longer than Ifm, since it breaks 
into pieces, as illustrated in Fig. [T2](b). The qq pairs produced from vacuum via the 
Schwinger mechanism completely screen the field of the end-point colour charges due 
to the linearity of the string potential [38]. The Schwinger phenomenon and existence 
of light quarks are the main reasons for not observing free quarks and gluons (colour 
screening) . 



5. Hadronization of colour charges 

Thus a colour charge is always accompanied by an anti-charge neutralizing its colour. 
The colour field in between forms a tube/string which is a very unstable construction, 
qq pairs pop up via tunnelling from vacuum, as shown in Fig. [T2r b). and the string 
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q 




Figure 13. Time-coordinate development of string fragmentation in the center of 
mass (left) and target rest (right) frames. 



is never much longer than 1 fm. The probabihty of such a string breaking over time 
interval T is given by, 



L J 

where L{t) is the time-dependent length of the string, and the probability density for 
the creation of a qq pair per unit time per unit length is given by the Schwinger formula 



The string length L{t) is getting shorter after each break, thus delaying the next 
pair production. Therefore, hadron momenta rise in geometric progression, i.e. the 
rapidity distribution of produced hadrons is constant. Notice that such a plateau in 
rapidity has been predicted by models with multiperipheral dynamics and for gluon 
radiation in perturbative QCD [39]. This process is illustrated on a time-coordinate 
plot in the cm. frame of the initial qq pair (e.g. e"'"e~ annihilation), and in the target 
rest frame (e.g. in DIS) in Fig.[T3]on the left and right correspondingly. Since both ends 
of the string are moving in the same direction and with the same velocity (the speed 
of light), the length of the string is independent of time. Its maximal possible value is 
Lmax = rnql However, after each break of the string it becomes about twice as short, 
as illustrated in Fig. [131 

Notice that the leading quark loses energy at a constant rate, dEq/dz = —n, though 
the whole hadronization process, until the creation of a hadron that includes this quark. 
It is interesting to notice that in perturbative QCD the leading quark loses energy for 
gluon radiation also at a constant rate, dEq/dz = —{2as/37i)Q'^ [40] . 




(35) 





(36) 
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Figure 14. Regge trajectories for mesons and soft Pomeron. 



6. Regge phenomenology 

The theory of Regge poles is a quite dormant topic. It does not seem to be taught very 
much anymore. In addition there is often found an attitude that the subject is obsolete, 
because it is identified so strongly with the pre-quark, pre-parton era of the S-matrix, 
dispersion-relations approach to strong interactions. This point of view is just plain 
wrong. The Chew, Frautschi, Regge, et. al., description of high energy behaviour in 
terms of singularities in the complex angular momentum plane is completely general. 
And the basic technique of Watson- Sommerf eld transform should be a standard part of 
the training in theoretical particle physics. 



6.1. Poles in angular momentum plane 

The energy dependence of the amphtude is governed by poles (or cuts) in the complex 
angular momentum plane [12] . 



where we sum over different Regge poles r, and hr{t) is a phenomenological residue 
function which is not given by the theory, but is fitted to data. It depends on t, but not 
energy, and correlates with the choice of the parameter sq. 

The phase factor ^^if) depends on the Regge pole signature a = (—1)'^, where J 
are spins (even or odd) of mesons lying on the trajectory. 



The energy dependent factor (s/sq)"*-*'' is controlled by the Regge trajectory a{t) 
which is nearly straight, a{t) = a{0) + as is demonstrated in on the Chew- Frautschi 



James Bj or ken [IT] 




(37) 




(38) 
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Figure 15. Dependence of the total cross-section (left) and elastic slope (right) on 
center of mass energy. Open and closed points correspond to pp and pp collisions, 
respectively. The curves show calculations based on the Pomeron model of 03] . Data 
are from Ref . [H] . 



plot in Fig. [TH This is the miracle of Regge theory: the linear Regge trajectories bridge 
the low-energy physics of resonances (t = > 0) with high -energy scattering (t < 0). 

High energies are dominated by Reggeons with highest trajectories arit), the 
Pomeranchuk pole (Pomeron), 

o;jp(0) ^ 1.1; 

a'jp ^ 0.25 GeV"^ , (39) 
and leading Reggeons, 

«/(0) ^ a^(0) ^ «p(0) ^ aaM ^ 0.5 ; 

a'j^ = 0.9 GeV-^ . (40) 

The first important prediction of the Regge pole theory was shrinkage of the elastic 
slope with energy. The slope parameter controls the t-dependence of the elastic cross- 
section, daei/dt oc e^*. According to (j37l) the slope parameter B rises with energy 

as, 

5(s) =5o + 2a^ln(s/so), (41) 

where Bq is a phenomenological parameter. 

The Pomeron parameters Eq. (139|) were extracted from data on elastic scattering. 
The Pomeron intercept ajp(O) comes from data on total hadronic (mostly pp and pp) 
cross-section fitted with the energy dependence (|37ll . while the parameter a'p is related 
to the elastic slope Eq. fHTl) . Corresponding data are shown in Fig. [151 
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Figure 16. The cross-section of single diffraction, a + b 
excitation channels at fixed effective mass Mx- 



P = 

bbb'b bbbb 

X + b summed over all 



6.2. Triple Regge phenomenology 

The cross-section for the inclusive process, a + b —>■ X + c can be also expressed in 
terms of the Regge approach. Here we focus on the most interesting case of diffractive 
excitation, c = b, via Pomeron exchange. To sum up all final-state excitations X, one 
can apply the unitarity relation to the pomeron-hadron {IP — a) amplitude as shown in 
Fig. [ini Provided that the effective mass of the excitation is large (but not too much). 
So <C Mjr -c s, one can describe the Pomeron-hadron elastic amplitude via Pomeron 
or secondary Reggeon exchanges in the t-channel. Then one arrives at the triple-Regge 
graph. Fig. [161 which corresponds to the cross- section, 

/ e \ Or(0) — 1 

E G^iPrW(l-XF)"'-(°)^'°"(*M-) , (42) 



dxp dt 



r=lP,]R ^-^U/ 

where xp is the Feynman variable for the recoil particle b defined in the center of mass, 
XF = 2pl/^s ^ l-Ml/s. 

Equation fj42l) contains new phenomenological functions, effective triple-Regge 
vertices, Gpjppit) and Cppjait). The diffractive cross-section can also be expressed 
in terms of the Pomeron-hadron total cross-section cr^"(s' = M^). Most interesting 
is the asymptotic (s' 
coupling, 

G^t) 



Ml > s) of this cross-section related to the triple-Pomeron 



(43) 



Here N]pbb[t) is the Pomeron-hadron vertex known from bb elastic scattering. Thus one 
can extract from data on single diffraction the Pomeron-hadron total cross-section, 
[l5] , which carries unique information about the properties of the Pomeron. The results 
shown in Fig. [T7| demonstrate an amazingly small cross-section, less than 2 mb. 

This is at least an order of magnitude less than one could expect. Indeed, the 
Pomeron as a gluonic object should interact more strongly than a meson, i.e., the 
Pomeron-proton cross-section could be about twice as big as the pion-proton one. Such 
a weak interaction of the Pomeron is probably the strongest evidence for the location 
of the glue in hadrons within small spots [16]. 
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Figure 17. The Pomeron-proton total cross-section extracted from single-diffraction 
data, pp — > pX^ as function of the invariant mass Mx which is the center-of-mass 
energy in IPp collision. Experimental data are from [45j . 



6.3. Building the Pomeron 

It has been a natural and simple assumption made in the early years of the Regge theory 
that the Pomeron is a Regge pole with a linear trajectory and the intercept ajp(t) = 1. 
Nowadays, however, we have a multi-choice answer, and it is still debated whether the 
Pomeron is: 

• a Regge pole (probably not, since ajp(O) varies with in DIS); 

• the DGLAP Pomeron [m HE] , which corresponds to a specific ordering for radiated 
gluons in the ladder graph in Fig. [271 vl+i < vl (see Section 7.2); 

• the BFKL Pomeron [19] which does not have ordering in transverse momenta of 
radiated gluons, but has no evolution with Q"^ either [50] (see Section 8); 

• something else? 

Gluons seem to be the most suitable building material: already the Born graph 
provides ajp(O) = 1. The higher order corrections are expected to pull the intercept 
above one. These corrections are dominated by ladder type graphs shown in Fig. [181 A 
ladder is a shadow of gluon bremsstrahlung according to the unitarity relation. Fig. [391 

The leading-log approximation (LLA) corresponds to keeping those terms only, 
where each coupling has a big factor ln(s). For fixed coupling, the BFKL result is 
not a Regge pole, but a cut with an intercept 




Figure 18. Perturbative Pomeron represented by the Born term, two-gluon exchange, 
and higher order terms having a form of gluonic latter graphs. 
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This result will be derived in Section 8, see Eq. (15^ . Unfortunately, the next-to- leading- 
log corrections (extra powers of Ug) to the intercept are of the same order |51j . 

«ip(0)-l= ^ ln2, (45) 

7r(l — Q.bas) 

and it may be even negative for most reasonable values of a^. 

However, it does not look reasonable to describe a soft Pomeron, controlling soft 
hadronic interaction at high energies perturbatively. Similar latter graphs, but built of 
light hadrons, e.g. of pions and a mesons as depicted in Fig. [19], well describe many 
features of soft hadronic collisions [52] . One can adjust the poorly known a-pion coupling 
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Figure 19. Ladder graphs built of pions and cr-mesons as a model for the soft Pomeron. 



to reproduce the Pomeron intercept. However, its closeness to one, which is very natural 
in QCD, looks as an accidental coincidence in this model. 



6.4- Duality 

Reggeons correspond to the exchange of valence quarks. The descriptions of meson- 
meson scattering amplitude in terms of interacting qq pairs in the t channel (Reggeons), 
or in the s channel (resonances) are dual [12], as is illustrated in Fig. [201 




Figure 20. The amplitude of meson-meson interaction via quark exchanges. Dashed 
lines show intermediate interactions within a quark- antiquark pair in i-channel (left), 
or in s-channel (right). 

No s channel resonances is possible in pp and K~^p elastic amplitudes. However, t 
channel Reggeons are present. To comply with duality the Reggeons must cancel each 
other in the imaginary part of the amplitude. For this reason, pairs of leading Reggeons 
must be exchange-degenerate, / with a;, and p with 02, i.e., their Regge trajectories and 
residue functions must be identical, differing only in the signature factors (phases) [12] . 
Data depicted on the Chew-Frautschi plot in Fig. [TH indeed confirm this expectation. 
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Figure 21. Electron- proton deep inelastic scattering 



The sums, f + uj and 02 + p must be real for pp and K^p, but imaginary for pp and 
K~p. Data at low energies dominated by Reggeons nicely confirm this. For the same 
reason spin effects are much stronger in pp and K^p^ than in pp and K~p. 

7. Deep inelastic scattering 

The basic idea of electron-proton Deep-inelastic scattering (DIS) is to use a lepton probe 
to study a hadron. A lepton of momentum k acquires momentum kl by exchanging a 
virtual photon of momentum q with the proton (with a spin 1/2) of momentum P and 
mass TTLj^ . After the collision, the rest of the energy is transferred to the unobserved final 



state X with mass M^. We ignore the lepton mass compared to the scale = — g^. 
The kinematics of DIS is characterized by a few Lorentz-invariant variables, see Fig. HTJ 
P -q 



V 

W'^iP + qf, 

s ={P + kf. (46) 
We define two other commonly used variables, namely the Bjorken variable 

= = ~ (A7) 

^ 2P-q 2mNy ^ + W^' ^ ^ 

where = —q"^, and the relative energy loss of the lepton 

p.q Q^ + W' 

y = p-k'^^^- ^^'^ 

The kinematic limits are < and x > Q'^/W'^ which leads to < x < 1. The value 
of X = 1 is reached when the proton is scattered elastically. The differential cross-section 
for inclusive scattering has the following form 

da{eP e'X) = Ysj2^> E(2^)'^'(^ + k - px - k')\A\\ (49) 

where one sums over all final hadronic states. The matrix element squared after summing 
over polarization of the virtual photon becomes 

1^1' = ^^(P|J^nO)|X)(X|J^(0)|P)V, (50) 
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where the leptonic tensor is 

L'^'' = {u{y)-i^'u{k)u{k)Yu{kt)), 

= lik^ky + k^k^ Q^iuk • k) . (51) 

In the above expression, we ignored the electron mass. We define the hadronic tensor 
as 



X 



= J d^xe'''^P\Jf'{x)r{0)\P), (53) 

where the second equation is obtained by using completeness for states of X. The 
hadronic tensor is directly related to the imaginary part of the forward Compton 
scattering amplitude via the optical theorem. We should stress that the hadronic 
tensor cannot be computed by perturbative QCD. However, one can write down the 
most general tensor from available momentum vectors P'^, q'^ and from g^'^ by using 
the transversality of the electromagnetic current q^W^" = q^W^" = , and parity and 
time-reversal symmetry W'^'^ = W^^: 

where Fi^2(a;,Q^) are the so-called structure functions. Making use of the above 
expression, the DIS cross-section reads 

- ^U-y- + ,^F,(x, g^)| . (55) 



In the next section, we consider whether Fi and F2 can be independent. The DIS can 
be viewed as 7*p scattering. The cross-section for a virtual photon with helicity A can 
be defined as 

a, = ^^e,{X)e:{X)lmT^\ (56) 

where the forward virtual Compton scattering amplitude T^'^ is related to the hadronic 
tensor, 

W^"" = —ImT^''' = lim—m^iqo + ie) - T^'^iqo - ie)), (57) 

2tI e^O 4:711 

with 

T^"" = ij d^xe'''^P\T{J>'{x)r{0))\P). (58) 

Here, T is the time ordering operator. Note that the relation between the hadronic 
tensor and the Compton scattering amplitude is a manifestation of the optical theorem. 
Using Eqs. ( |56]l58l) . one obtains the corresponding cross-section for transverse and 
longitudinal photons. 



Q^{1 - X) 



Q2(l-X) 



1 + 91\f2{x,Q')-2xF,{x,Q') 



(59) 
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It is sometimes common to define linear combinations of the structure functions 

Fl{x, Q2) = F2{x, Q^) - 2xF,{x, Q^), (60) 
Ft(x,Q2) = 2xFi(x,g2). (61) 

The usefulness of the above definition is that 7*p scattering for transverse and 
longitudinal photons can be defined in terms of -Fl^t- 

7.1. Bjorken scaling and parton model 

In the late 1960's, experimental result from SLAG [53j surprisingly showed that the 
structure function F2{x,Q'^) is nearly independent of at a fixed x. In Fig. [251 we 
show the measured values of F2{x,Q^) for various as a function of x. It is obvious 
that all data points seems to lie on a single curve which show that within error bars F2 
is independent of Q^. This phenomenon is called Bjorken scaling \5^ 

An intuitive explanation of this phenomenon was given by Feynman [56] , within the 
parton model. In the parton model, one assumes that the proton is made of pointlike 
charged constituents, called partons which interact incoherently. Then the total 7*p 
cross-section can be written as an incoherent sum of photon-parton cross-sections. We 
work in the Breit frame in which proton and virtual photon are moving collinearly 
and the virtual photon does not carry the energy, but only momentum. Assume that 
the scattering is elastic and the parton of type q carries a fraction rj of the proton's 
momentum, see Fig. [231 For massless partons, we have 

(g + r]Pf = 2r]P.q - = 0, (62) 

which leads to 

T] = X. (63) 

This implies that in the Breit frame, the Bjorken x is the momentum fraction of the 
proton carried away by the struck quark. 

In order to calculate the 7*p scattering cross-section in the parton model, one should 
first calculate eg — >■ eq cross-section, which can be obtained from those for e"'"e~ — >■ qq by 
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" xP+q 



(l-x)P 

Figure 23. Parton picture of DIS and Bjorken x. 

crossing symmetry. Equivalently, one may first calculate the cross-section for transverse 
and longitudinal photons scattering off spin-1/2 parton, see Fig. [23l 

= 0. (64) 

For massless quarks, the longitudinal cross-section has to be zero because of helicity 
conservation. Now, by comparing Eqs. fl64|59l) . one can define the structure function in 
parton language, namely by introducing the density qf{x) of quarks of flavor / inside 
proton, 

F2ix)=x Z}{qj{x) + qj{x)), (65) 

f=u4,.. 

Fl = 0. (66) 

At this order, the structure functions depend only on x and not on Q^. When the 
longitudinal structure function vanishes, one obtains the Callan-Gross relation [58] , 

F2 - 2xFi = 0. (67) 

This equation is approximately confirmed by experiment and proved that partons are 
fermion with spin one half. It is common to define valence quarks Uv,dv, ■■, as 

dv =qd- Qd, (68) 

where g„, qd are called sea anti-quarks. It is also possible to measure DIS on the neutron 
and extract the neutron structure function. Assuming strong isospin symmetry, we have 
the following relations between the parton distribution functions in proton and neutron: 

qu/n = qd/p, qu/n = qd/p, qd/n = qu/p, qs/n = qs/p, qc/n = qc/p, 

(69) 

and so on. The convention is always to refer to the proton PDFs and drop the subscript 
/p or /n. 

At this stage, one may wonder why the proton form factor F{Q) falls steep with 
Q, while the structure function does not. The answer is that the form factor is the 
probability for the proton to survive intact a kick of strength Q. The stronger the kick, 
the less survival probability. However, in the case of inclusive DIS, all final states are 
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Figure 24. Comparison of the measured proton structure function F2 with QCD fits 

m- 

allowed, so the total probability saturates and is independent of Q. A similar situation 
is known to happen in hadronic collisions, where the t-slope of single diffraction is half 
that for elastic pp, because of the disappearance of one of the proton form factors. 

7.2. Scaling violation and DGLAP evolution equation 

In the previous section we showed that at leading order the partonic sub-process of DIS 
eg — i> eg is independent which leads to Bjorken scaling. This could be correct if the 
number of partons were constant. However, they are not classical particles but quantum 
fluctuations. A photon of virtuality Q can resolve partons with transverse momentum 
< Q but is blind to harder fluctuations. Increasing Q, one can see more partons in 
the proton. Correspondingly, the parton distribution slowly changes with Q, shifting to 
smaller x due to momentum conversation, i.e. it is expected to rise with Q at small x, 
but to fall at large x, see Fig. [211 The dependence of the structure function can be 
described by the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) equations [17]. 

One of the major successes of QCD has been the prediction of the pattern of Bjorken 
scaling violation as shown in Fig. [211 We will explain that the presence of gluon radiation 
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Figure 25. Four type of diagrams corresponding to DGLAP splitting functions of 
QCD. 



controls the behaviour of Bjorken scahng violation. At higher order in one should 
also include gluon radiation eq eqg. As in the previous section, one obtains at next 
order 




-In 




(70) 



where g{x\) denotes the gluon density of the proton. The origin of the ln((5^//i^) is easy 
to understand. The struck quark acquires transverse momentum pt with probability 
tts^-i^. On the order hand, partons with p\ > are suppressed. Now, integrating over 
all phase space, pt produces the logarithmic term aslnQ^//!^. The parameter fi was 
introduced as a cutoff regulator. The divergence when /i corresponds to a case that 
the outgoing gluon becomes exactly collinear with the incoming quark. This means that 
the internal quark line becomes on-shell leading to the logarithmic divergence. This is 
called collinear divergence. The function Pfg is quark-quark splitting function [171 ES] 

..W^i(i^). (n) 

The splitting function Pfg shows the probability for a quark to turn into a quark and a 
gluon. This function is independent of the regularization and is universal. 

Note that the structure function equation (ITOll is now obviously dependent, 
violating Bjorken scaling. We also introduced an ad hoc parameter fi, called the 
factorization scale, which separates the underlying physics into two parts: all physics at 
scale below /i contained in the parton distribution, and all calculable physics at scales 
above n are part of the partonic scattering cross-section. It is important to mention 
that although we have obtained Eq. (170|) at the higher order, the leading logarithmic 
behaviour is universal and this factorization formula is valid at any order of as- Since 
fi is not a physical quantity, observables should not depend on it. Therefore, 

rfF2(x, Q^/i) 
dfi 

and 



0, (72) 



^2 dqf{x,Q^) as dxi fx., , , 

Q 17T2 = ^ / —PfG — 9{xi, Q )• (73) 



dQ"^ 271 J X Xi \Xi 
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This is one of the DGLAP equations [17] which describes the evolution of the quark 
density. In the same fashion one can obtain the DGLAP equations for gluon density 
giyXi^Q'^). Altogether one obtains Nf + 1 coupled equations (ignoring anti-quark for 
simplicity) describing the evolution of the singlet parton densities (?/(a;i,(5^) and 

( ^/(-' Q') \ a. } dx, ( (t) (^) ^ ( ^/(-- 



^ dQ^ 



[ g{x,Q') I 27r7 ^ p^^^^^ j y ^^^^^q 



(74) 



where the splitting function Pff,Pfg,Pgf and Pgg can be computed from pQCD order 
by order. The analytic calculation of these splitting functions to next-to-next-to-leading 
order has been carried out in Ref [SH] . In Fig. [25] we show the lowest typical diagrams 
corresponding to the various splitting functions. 

The steep rise of F2 at small x in Fig. [21] can be simply seen from the double log 
DGLAP equation, 

d'^xg{x,Q^) N^as , /vp-n 

mn(i/x)mng^ = ^"^^"'^^- ^''^ 

For a fixed coupling constant a^, the solution can be approximated by 



xgix, g2) oc exp |^2^ ^ ln(l/x) In(QVQo) j • (76) 

This equation clearly indicates that at small x and high Q'^ the gluon density rises. 

The parton distribution function (PDF) cannot be calculated from the first 
principles. However, their scale evolution can be perturbatively computed via DGLAP 
equations. We therefore, calculate the fj,"^ dependence of the PDFs. In this way knowing 
the value of PDF at a given scale by fitting data is sufficient to obtain information about 
PDFs at all scales via DGLAP evolution equations. The DGLAP equation is a special 
kind of a renormalization equation. It is obvious from the In Q"^/ /i^ term that one should 
not chose /i too far from since the log term will become large enough to compensate 
the smallness of as and perturbative computation will become questionable. 

The typical strategy for extracting PDFs from DIS data is first to introduce ad hoc 
PDFs at some scale and then to evolve them with DGLAP to other scales, and finally 
compare F2 at higher values of with data and adjust the starting PDFs. Having 
good data with high statistics one can single out PDFs for different parton species. 
Such parametrizations are provided by three collaborations: GRV [62], MRST [63] and 
CTEQ [6l] in leading and next to leading orders. In Fig. [26] we show typical PDFs at 
= 10 GeV^ as a function of x. 

Having results for the PDFs, one can check how much of proton's momentum is 
carried by quarks and antiquarks. The data shows that, 

dxF2 ^ 0.5. (77) 



This result is quite significant, since it shows that only half of the total momentum 
is carried by all quarks and antiquarks in the proton. Another half of the proton 
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Figure 26. PDFs xf{x,^P) at 20,10"* GcV^ as a function of x. It is clearly 

seen that the gluons and the valence quarks are more important at small and large x, 
respectively. The curves are obtained from the NNLO global analysis [60]. The figure 
is taken from Ref. [ST]. 

momentum is carried by partons which do not interact with the photon, apparently 
gluons. 

7.3. Factorization theorem 

A cross-section of any hadronic reaction with a hard scale generally gets contribution 
from short- and long-distance interactions, and is hence not computable directly in 
perturbation theory for QCD. Factorization theorems [65] allow one to derive predictions 
for the hadronic cross-sections by writing the cross-section as a convolution product 
of factors, namely an infrared finite part for the short distance which is calculable 
in perturbative theory, with a nonperturbative function which is universal to many 
different processes, but non-calculable at the perturbative level. The applications and 
predictability of perturbative QCD rely on the factorization theorem. 

There has been tremendous effort to examine factorization theorems for various 
processes; for a review see Ref. [65]. For inclusive processes, it has been shown that 
the factorization theorem holds if (1) all Lorentz invariants defining the process are 
large and comparable, except for particle masses, and (2) one counts all final states 
that include the specified outgoing particles or jets, namely in processes as hadron 
A + hadron B hadron C + X, the X denotes anything else, in addition to the 
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Figure 27. The QCD improved parton model, Xi denote the longitudinal momentum 
fraction of the partons with respect to the target. 



specified hadron C. 

For example, in DIS, the factorization tlieorem for the structure functions has the 
following form. 



which is valid in the Bjorken limit in which Q gets large with x fixed. The sum is over 
all species of partons, namely gluon, quarks and antiquarks of different flavours. The 
function fa/H denotes the PDF of parton of type a in hadron H. The hard process- 
dependent factor Cia is ultraviolet dominated, that is, it receives important contributions 
only from momenta of order Q. This ensures that one can perturbatively calculate Cia 
in power of as{Q) [see Eq. (170|) ]. Notice that the factor Cia depends only on the parton 
type a, and not directly on our choice of hadron A. The parameter /i in Eq. (ffSll defines 
the limit between the short- distance dynamics. The ability to calculate the Cia leads to 
great predictive power for factorization theorems. For instance, if we measure F2{x, Q) 
for a particular hadron A, Eq. ( ITHl) will enable us to determine the PDFs fa/ a- Then we 
predict Fi{x, Q) for the same hadron A, in terms of the same fa/ a and calculable Cia- 

8. BFKL formalism 

We recall that the DGLAP equations take into account all the contributions proportional 

to 



which arises from ladder type diagrams with strong ordering in the transverse momenta, 
see Fig. [271 i-e., 



For processes where x is small, but is not sufficiently large to make the double 
logarithmic approximation valid, the Balitsky-Fadin-Kuraev-Kuraev-Lipatov (BFKL) 




(78) 




(79) 



Pti » Pt2 » ■■■■ 



(80) 
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equation [IS] has been proposed. In this scheme, the gluonic branching in the ladder 
diagrams has ordering in longitudinal momentum (see Fig. [27|) 

Xi >> X2 » .... (81) 

This resums, in the so-called ln(l/a;) approximation, the terms 

[a.(Q2)lnf-V- (82) 



At the same time there is no ordering in transverse momentum, one may have 

Pti ~Pt2 ~ •••• ~PTn- (83) 

In the high energy limit, the scattering processes are dominated by partonic processes 
with gluon exchange in the t-channel. The BFKL equation accounts for resummation 
of multiple gluon radiation when s » t. 

The BFKL equation is more conveniently written in terms of the unintegrated gluon 
density (j){x, fcf.) which relates to the gluon density g{x, Q"^) introduced in the previous 
section by 



?2 

xg{x,Q'')=l dkl<p{x,kl). (84) 
Jo 

The unintegrated gluon distribution gives the probability of finding a gluon in the hadron 
with longitudinal momentum fraction x and transverse momentum kx- Note that there is 
no unique definition for the unintegrated gluon density in terms of gluon density 
For comparison of various parametrizations for the unintegrated gluon distribution in 
different schemes, see Ref. [66] . 

At leading order in ln(l/x) the BFKL equation can then be written in the following 
simple form: 

(90(x, kl) _ Ncas f dpi ( ^^^^ _ kl(f){x, k^) \ ^^^^ 



dhi{l/x) Tc"^ J {kr — PtY \ ' Pr ^ i^T — PtY / 

Equation ( |85l) is illustrated in Fig. [281 The first term in Fig. [28] corresponds to two 
gluon exchange, the initial condition for Eq. fl85l) . The first and second terms on the 
right hand side of Eq. ( l85i) , correspond to the second (real) and third (virtual) terms in 
Fig. [281 respectively. Iterating the BFKL kernel leads to the ladder diagrams shown in 
Fig. [IHI Note that the gluon propagators and vertices in Figs. [T8ll28l are not the usual 
QCD vertices and propagator. The vertices are effective Lipatov vertices and propagator 
are the so-called reggeized gluon propagators generated by iterating the BFKL kernel 
[19] . We refer the interested readers to Refs. [l9l [68] [69] for derivation of the BFKL 
equation. For a recent review of the subject see Ref. [70] . 

It is rather straightforward to see that for a fixed coupling a^, the solution for the 
unintegrated gluon density up to a constant is 

0(a;, 4) oc (l/x)"^-\ ajp-l = ^^^ln2. (86) 

TT 

However, it is important to notice that, based on the BFKL equation, number of gluons 
rises with 1/x forever. This strongly indicates that some physics must be missing here. 
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Figure 28. Schematic representation of BFKL evolution for unintegrated gluon 
distribution. The dashed vertical denotes the cut. 



That is because in QCD the gluon fields cannot be stronger than ~ l/gf at very 
small coupling g. Therefore, when the gluon field reaches a density with 
F 1 

^^--.^ (87) 



we expect some new physics to be at work in order to slow down the rise of gluon density. 
We shall postpone to elaborate more on this problem in the next section. 

The total cross-section of quarkonium-quarkonium scattering, in the lowest order in 
the coupling a^, in the simplest model of two-gluon exchange between two quarkonium 
is energy independent 

a ~ s°, (88) 

where s denote the center of energy of the system. However, experimental data indicates 
that hadronic cross-sections increase as power of s, see Fig. [15] 

a ~ s^, (89) 

where A is called the Pomeron intercept. Finding an explanation for the experimental 
value of A has been one of the remaining challenges of QCD. The DGLAP equations has 
been successful in describing the DIS data, but it cannot reproduce the energy growth of 
hadronic cross-sections. One of the interesting features of the BFKL formalism is that it 
naturally leads to an energy-dependent cross-section. One can show that cross-sections 
mediated by the BFKL exchange grows as a power of energy 

a ~ s"^-\ (90) 

where ajp is given in Eq. ( l86l) . Unfortunately the value of Oip — 1~0. 8 is higher than 
experimental value 0.2 — 0.3 observed in DIS experiments. Before going to higher-order 
corrections to the BFKL kernel, it is important to notice that there is already a serious 
problem at lower order. That is due to the fact that the power energy growth of the 
total cross-section fl90|) violates the Froissart unitarity bound [71] which put a limit on 
the growth rate of total cross-sections with energy s at asymptotically high energies 

a < const In^ s. (91) 

This indicates that the BFKL kernel should be modified in order to restore the unitarity 
at high energy. 
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9. The GLR-MQ evolution equation and saturation 

In the previous section we pointed out that, based on the BFKL formalism, the number 
of gluons rises sharply at small x or high energy. At the same time, the transverse 
sizes of the gluons tt ~ 1/pt, can be similar, see Eq. ( l83l) . This means that at high 
energy, a hadron produces many gluons with a similar size. As the energy increase, 
more gluons are produced and eventually they start overlapping in transverse plane. 
The crucial assumption behind both DGLAP and BFKL evolution equations is that 
parton densities inside a hadron are small enough, so that the only important partonic 
sub-process is splitting. However, at very low value of x, the gluon density may become 
so large that gluons start overlapping and the gluon recombinations process becomes 
important. This phenomenon is generally known as parton saturation, and it should 
limit the growth of the gluon density generated by splitting. 

Gribov, Levin and Ryskin (GLR) [72] proposed that at high density of gluon 
fields when nonlinear effects become important, there should be an energy region where 
the gluon recombination becomes important. In the GLR scheme this recombination 
is described through a modification of the linear BFKL equation with a quadratic 
correction which gives rise to effective ladder merger vertices which are the triple 
Pomeron ones, 

dct>{x, k\) N^a. f dpi ( _ kU{x,kl) \ _ ^^^^^^ ^^2) 



91n(l/x) 7r2 J {^kT-pxYy^ " pl + {kT-pTY) St 

where St = TfR^ defined the geometrical cross-sectional area of a hadron or a nucleus 
along the beam axis. In comparison with the BFKL equation (185!) . only the last term 
is new. 

Later, it was proved by Mueller and Qiu [73] that the ansatz Eq. fl92l) can be 
derived in the double leading logarithmic approximation (DLLA) with a resummation 
of the type (a, ln(QVA2) ln(l/x))". Muller and Qiu [73] showed that in the DLLA 
approximation, including diagrams with two fusing DGLAP ladders, one arrives at the 
following nonlinear equation for gluon density: 

The above equation can be converted into Eq. (!92l) via the definition of an unintegrated 
gluon density Eq. Eq. (jHSD is known as the GLR-MQ equation. 

Notice that in the DLLA approximation both the BFKL and DGLAP equations are 
identical, since the resummations are the same. This can be already seen in Eq. (1751) . 
where the merging ladders were ignored. That equation is identical to the first term of 
Eq. (I93D. 

One of the remarkable properties of the GLR-MQ equation fl93|) is that it introduces 
a scale at which the non-linear effects become relevant. This may occur when the 
linear and quadratic terms in Eq. ( 193|) becomes equal: 

s^.. — ■ ^^^^ 
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A quantitative condition for gluon saturation can be obtained by comparing the 
gluon recombination cross-section a ~ OLsjQ^ with the surface density of gluons 
p ~ xg{x^ Q'^)/'kR^. Saturation takes place when crp ~ 1 which leads to Eq. ( 19^ . 

The saturation scale separates the linear (governed by DGLAP or BFKL 
equations) and non-linear evolution of QCD. The DGLAP, the BFKL and saturated 
regimes are sketched in Fig. [291 At low energy, colour screening is due to confinement 
with typical colour screening distance ^qq£)) ^-^id thus non-perturbative. At high 
energy (or small x), partons are much more densely packed, and colour neutralization 
occurs in fact over distances of the order Qg « A~^. This means that small x physics 
seems to be universal, and all hadrons and nuclei should behave in the same way at very 
high energy. 

The basic physics of saturation is to introduce higher twist terms [721 ES] in the 
factorization formula like Eq. (!70l) . This is difficult to implement. During last decade 
there has been some progress along these lines and some models has been proposed 
|70] . For example, the description of this non-linear evolution has been given in the 
so-called Colour Glass Condensate [71] scheme in terms of a classical field theory of 
dynamical gluon fields coupled to static stochastic sources. The evolution of multi- 
parton correlators with energy is described by the JIMWLK renormalization group 
equations [75]. At large Nc and large nuclei, one recovers the Balitsky-Kovchegov (BK) 
equation [76] for forward colour dipole cross-section. The recombination effect is taken 
into account by the non- linear term of the BK equation [76l [7^ . 
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Figure 30. Photon virtual dissociation to a qq pair with transverse separation 
and fractional the light-cone momenta a and 1 — a. 



10. The colour dipole approach and low-x DIS 

The parton model description is not Lorentz invariant, only observables have to be 
Lorentz invariant. One cannot even say where a sea parton has originated, who is the 
owner, the beam or the target, see Fig. [301 In the domain of small x, sea quarks and 
gluons dominate, and the rest frame of the proton is more convenient. In this frame, the 
photon can convert into a quark-antiquark pair which then develops a parton cloud, see 
Fig. [301 Photons can also hit a quark inside the target without being split to gg-pair. 
However, in target rest frame this process is strongly suppressed. Therefore the former 
contribution is the dominant one. The lifetime tc of such gg-pair fluctuation can be 
estimated via the uncertainty relation tc ~ 2m^^ where niN = ^ GeV is the mass of 
a nucleon. The smaller the Bjorken x, the larger the coherence time. For the lowest 
value of X accessible at HERA the coherence time in the proton rest frame is about 10^ 
fm. Therefore the coherence time or lifetime of such a pair creation tc can be larger 
than nuclear radius at low x and pairs can experience multiple scattering within the 
coherence length. This is very important point for understanding the phenomenon of 
nuclear shadowing. The total 7*— p cross-section, or the forward amplitude, is described 
as the interaction of a qq fluctuation of the photon with the target, as is shown in Fig.[30| 
The cross-section for the transverse and longitudinal photons is given by the 
factorized formula, [251 [78] , 







all = I da I d'vT ^lf{a,TT)' a,q{TT), (95) 



2 



where rj- is the distance between the quark and antiquark in the transverse plane and a 
is the fraction of the photon energy carried by the quark, see Figs. [301 and [HH The cross- 
section for scattering a gg-dipole off the proton is denoted by agq^rx). The light-cone 
(LC) distribution functions \&^^^(a,rT) for the transition 7* qq can be calculated in 
perturbation theory and read to first order in the QED coupling constant ac„ 



9 AT f 

^?.("' ^4 = 5: - 2a(l - a)] e-Kl{err) + mjK^(6r^)} , (96) 
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Figure 31. The virtual photon interacts via its hadronic fluctuations which are qq 
dipoles and more comphcated Fock states. The Pomeron exchange is illustrated as a 
perturbative ladder. 



12^ 



Nf 



/=i 



a, 



(97) 



where Ko,i are the modified Bessel functions of the second kind. Note that the above 
distribution functions are not normahzed and can be even divergent. That is why we 
avoided to call them wave functions. Although the transverse part of the distribution 
function is divergent at 0, the dipole cross-section vanishes in this limit as 

(^qqirr) oc rf,, so the result of Eq. ( 195|) remains finite. We have also introduced a 
parameter 

= a(l-a)Q^ + mJ, (98) 

where the parameter is quark mass. The mean transverse qq separation for a virtual 
photon is controlled by the Bessel functions, 

^^^^ = ail-a)Q^ + my ^^'^ 

Thus the separation is about as small as except at the end points a — > 0, 1. This 
implies that even a highly virtual photon can create a large qq fiuctuation although with 
a small probability. This is an important point for the aligned jet model [80]. Notice 
that ~ ^QCD P^^ys here the role of an infra-red cutoff. 

The incoming photon (or hadron) is not an eigenstate of the interaction, since it can 
be diffractively excited. Therefore one should switch to the eigenstate representation. 
The choice of the eigenstate basis depends on the underlying theory. It was first realized 
in Ref. [25] that the eigenstates of interaction in QCD are colourless dipoles. Such dipoles 
cannot be excited during the interaction time and can experience only elastic scattering. 
Indeed, high-energy dipoles have no definite mass, but only separation which cannot 
be altered during soft interaction. The eigenvalues of the total cross-section aqqirr) 
depend on r^, but may also depend on energy. 

At the level of two-gluon exchange (Born approximation), the dipole cross-section 
is independent of energy and related to the two-quark form factor of the proton via [25] 



11 



(rr) 



16a? 



[1 - {p\ exp(ipT ■ (n - r2))\p)] [1 - exp(ipT ■ ^t)] 



4 

Vt 



(100) 



Notice the colour screening factor [1 — exp(i|>r ■ ^t)] in Eq- fHOOP - which makes the 
dipole cross-section vanishes as r\ at 0. This is an important property of the 

dipole cross-section which is the cornerstone of the colour transparency phenomenon. 
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Figure 32. Photon-proton total cross-section as a function of r = Q^/Q^g for a; < 0.01. 
The data are from the HI [85], ZEUS [86], E665 [87] and NMC |88] collaborations. 



The energy dependence of the dipole cross-section is generated by higher order QCD 
corrections. For small distances — >■ 0, one can relate (Jqq{rT) to the phenomenological 
gluon density [81] 

cr,,-(x, Tt) = y4a,(Q2 ~ l/r^T)^9{x, ~ 1/4). (IQl) 

When the dipole cross-section is proportional to the gluon density of the target, only 
quarks generated from gluon splittings are taken into account in the cross-section 
Eq. ( |95l) . In other words, the valence quark contribution (or the reggeons in the Regge 
phenomenology) are neglected and therefore Eq. fl95|) is only applicable when sea quarks 
dominate, i.e. at low x. Having said that, the master Eq. ( l95i) is quite general and does 
not rely on the applicability of the pQCD. 

The dipole cross-section is theoretically difficult to predict, but several 
parametrizations have been proposed in the literature. For our purposes, here we 
consider two parametrizations, the saturation model of Golec-Biernat and Wiisthoff 
(GBW) [82] and the modified GBW coupled to DGLAP evolution (GBW-DGLAP) 



10.1. GBW model 

In the GBW model [H2] the dipole cross-section is parametrized as, 

Mx,r) = ao(l-e-i^'«^(^)), (102) 

where the parameters, fitted to DIS HERA data at small x, are given by cxq = 23.03 
mb, Qs{x) = IGeV x (x/xo)"^/^ where Xq = 3.04 x 10"^ and A = 0.288. This 
parametrization gives a quite good description of DIS data at x < 0.01. One of the 
interesting feature of the HERA data is a geometrical scaling [SI]; namely all available 
data for the inclusive virtual photon-proton cross-section for < 0.01 and various 
seems to scale as a function of r = Q^/Q^, see Fig. [321 This might indicate that the 
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Figure 33. Right: The dipole cross-section for x — 10^^, 10"'^, ...10~^ from left 
to right. The sohd and dotted hnes show results from the model with the DGLAP 
evolution Eq. (I103|) and the saturation model Eq. (|102p . respectively. Left: The 
effective slope \{Q^) from the parametrization F2 ~ x'^'^^ ^ as a function of Q^. 
The lines are the same as the right panel. The figure is taken from |83) . 



semi-hard scale which is also present in the saturation region (See Section 9), plays a 
role already at the kinematics of HERA. However, one should be aware that the DGLAP 
evolution describes the same data as well. So far it is not clear how much saturation is 
relevant to available DIS data. 

A salient feature of the model is that, for decreasing x, the dipole cross-section 
saturates for smaller dipole sizes, and that at small r, as perturbative QCD implies, the 
colour transparency phenomenon cr ~ r^, is at work. 

One of the shortcomings of the GBW model is that it does not match QCD evolution 
(DGLAP) at large values of Q^. This failure can be seen in the energy dependence of aj^t 
for > 20 GeV^, where the model predictions are below the data [82l[83]. Apparently, 
the simple rf. behaviour at small tt should be corrected. 

10.2. GBW coupled to DGLAP equation and dipole evolution 

A modification of the GWB dipole parametrization model, Eq. fll02p . was proposed in 
Ref. p]: 

7rVa,(/i2)x^(x,/i2)^^ ^^^^^ 




3(70 

where the scale /i^ is related to the dipole size by 

/^' = -^ + /io- (104) 
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Here the gluon density g{x,ii^) is evolved to the scale iJ? with the leading order (LO) 
DGLAP equation ([71]). Moreover, the quark contribution to the gluon density is 
neglected in the small-x limit, and therefore 



The initial gluon density is taken at the scale Q\ = IGeV^ in the form 

xg{x, n^) = AgX'^^ (1 - xf-^, (106) 

where the parameters C = 0.26, /ig = 0.52GeV^, Ag = 1.20 and Xg = 0.28 are fixed 
from a fit to DIS data for x < 0.01 and in a range of between 0.1 and 500 GeV^ 
[83] . We use the LO formula for the running coupling a.^, with three flavors and for 
Aqq£) = 0.2 GeV. The dipole size determines the evolution scale /i^ through Eq. fll04p . 
The evolution of the gluon density is performed numerically for every dipole size r. 
Therefore the DGLAP equation is now coupled to the master equation (jHSD- It is 
important to stress that the GBW-DGLAP model preserves the successes of the GBW 
model at low and its saturation property for large dipole sizes, while incorporating 
the evolution of the gluon density by modifying the small-r behaviour of the dipole size. 
Fig. [331 

To highlight the failure of GBW parametrization, in Fig. [33] we show the effective 
slope A((5^) from the parametrization F2 ~ x"'^*^'^^-' as a function of Q^. It is seen that 
the GBW-DGLAP parametrization is essential in order to describe the data. 

11. The Drell-Yan process and direct photons 

11.1. The partonic description 

As we already mentioned, the PDFs are universal. Therefore one can use the DIS data to 
extract PDFs and then make prediction for other hard processes. The most prominent 
example of hadron hadron collisions is the so-called Drell Yan (DY) process [89] , where 
lepton pairs are produced: 

/ll + /l2 /i+ + yU" + X, (107) 

where X can be any undetected particles. In the parton model, this process looks like 
a quark and an anti-quark from two hadrons annihilating into into a lepton pair, see 
Fig. [31 

The kinematic of the DY process can be conveniently defined via light-cone 

momentum fractions of the projectile (target) parton, xi {X2), 

m-q 2P,-q 
xi = , X2 = . (108) 

s s 

where Pi and P2 are the four momenta of hadron 1 and hadron 2, respectively and s 
denotes the square of the center-of-mass energy of the colliding hadrons is s = (Pi + P2)^- 
The Feynman variable xp is related to other kinematics variables as, 

Xf = ^Xi-X2, (109) 
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Figure 34. Partonic picture of the DY process in the leading order. Two hadrons 
colhde and a quark from one hadron annihilates with an antiquark from the other 
hadron into a timelike photon, which decays into a lepton pair. 



where p'f" is the longitudinal momentum of the dilepton in the hadron-hadron center- 
of-mass frame. Another relation is 

r = X1X2 = , (110) 

s 

Where = > denotes the mass of the spacelike photon and the transverse 
momentum of the virtual photon has been neglected. The partonic annihilation cross- 
section for Fig. [34] reads 

da _4-«e^™^|^(^^,^,_M^). (Ill) 



dM^ 3N^M^ 

The hadronic cross-section can be then written as the convolution of PDFs with the 
partonic cross-section, like in DIS, 

+ (1-2)};^. (112) 

where qf{xi) is the probability to find a quark of flavor / with light-cone momentum 
fraction xi in hadron a, and g/ is the analog for antiquarks. In the second line, we 
have plugged the partonic cross-section Eq. (II lip and performed one of the integrals. It 
is interesting to note that the right-hand side of Eq. (I112p depends only on r and not 
separately on and s. This scaling property is confirmed experimentally [90j ■ 

Some features of dilepton production cannot be understood in the lowest order 
picture. The cross-section given by Eq. (I112p is 2 — 3 times smaller than the measured 
value. This discrepancy is usually treated by introducing an ad hoc normalization factor, 
the so-called K factor. The K factor is approximately independent of M^. Another 
obvious problem is that the transverse momentum spectrum in the naive parton model 
cannot describe data. Phenomenologically, one can introduce a primordial momentum 
distribution of the quarks, but what is observed in experiments about 1 — 2 GeV is much 
larger than what one would expect from Fermi motion. 

These problems can be partially resolved by taking into account the next-order 
QCD corrections, shown in Fig. [351 Owing to the radiation of the gluon (the second 
row diagrams), the quark acquires a transverse momentum. In this way, the pQCD 
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Figure 35. Higher order QCD 
corrections to the DY process. 
The diagrams for virtual correc- 
tions, the annihilation process, 
and the Compton process are de- 
picted in the upper, middle, and 
last row, respectively. These 
higher order corrections account 
for most of the K factor and ex- 
plain data at large transverse mo- 
menta. The figure is taken from 
Ref. [S]. 



correction provides the missing mechanism for the production of lepton pairs with 
large transverse momentum px- However, the transverse momentum spectrum is not 
described well in this order, and obviously somethings is still missing. In particular, at 
low pt the pQCD result diverges. There have been attempts to overcome this problem 
by a resummation of soft gluons radiated from the quark and antiquark [02] • The last 
row in Fig. [35] displays the diagrams for the QCD Compton process, where a quark 
in one hadron picks up a gluon from the other hadron and radiates a photon. This 
mechanism is dominant at large px [93] . 



11.2. The colour dipole description 

Similar to DIS, the DY process can be viewed in the target rest frame where it looks like 
bremsstrahlung rather than parton annihilation, see Fig. [36l A quark or an antiquark 
from a projectile hadron radiates a virtual photon while hitting the target. This 
radiation can occur before and after the quark scatters off the target. The impact 
parameter representation of the cross-section for such a process can be written in the 
factorized form similar to DIS [HI [HSl 1^ . 

da{qp^qYp) f ,2 i,t,t,l^ ^|2 / n n 1 Q^ 

— = / d tt W i (a,rT) aqq{x,arT), (113) 

a ma J 

where a is the light-cone momentum fraction of the quark, carried away by the photon 
and rT the transverse separation between 7* and q. The dipole cross-section aqq{x, ar-r) 
with transverse separation ary is a universal quantity like PDFs and has already been 
introduced in the DIS section. 
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Figure 36. In the target rest frame, the DY process looks Hke bremsstrahlung. A 
projectile quark (or antiquark) scatters off the target and radiates a massive photon 
which subsequently decays into the lepton pair. The photon can also be radiated before 
the quark hits the target. Both diagrams are important. 



Where does the dipole cross-section come from if there is no dipole in diagrams 
given in Fig. [361 ^"^^ "^hy is the transverse dipole size ar-r? The dipole cross-section 
appears because the quark is displaced in the impact parameter plane after radiation of 
the photon. The antiquark enters after taking the complex conjugate of the amplitude. 
Therefore the dipole in equation flllSp is not a real gg-dipole. As in the real dipole in DIS 
where colour screening is provided by interactions with either the quark or the antiquark, 
in the case of radiation the two amplitudes for radiation prior or after the interaction 
screen each other, leading to cancellation of the infra-red divergences. Now, back to the 
second question, if tt is the transverse separation between the quark and the photon, 
and a is the fractional momentum of radiated photon. Then, the transverse separation 
between the photon and center of gravity is (1 — a)rT and the distance between the 
quark and center of gravity will be ar^. Before radiation, the center of gravity of the 
quark coincide with the incident quark, after radiation the relative distance between 
quark and the center of gravity is shifted to ar-p. Taking the complex conjugate of the 
amplitude, it looks as if the transverse size between q and q is arx which is the argument 
of the dipole cross-section. 

The wave function of the 7*g fluctuation in Eq. (11131) for transversely and 
longitudinally polarized photons reads, 

^^4(a,r^) = ^XfO^^'x.MwT). (114) 

Here Xij ^I's the spinors of the initial and final quarks and Kq{x) is the modified Bessel 
function. The operators O^'^ have the form, 

= imfO^e* ■ (n x a) + a e* ■ (a x V) - i(2 - a) e* ■ V , (115) 

= 2M(1 - a) , (116) 

where e is the polarization vector of the photon, is a unit vector along the projectile 
momentum, and V acts on r^- For radiation of prompt photons M = 0. A parameter 

r^2 = ^ct^ + (1 - a) , (117) 
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is the analog of the parameter e Eq. fl98p in DIS. 

In order to obtain the hadronic cross-section from tlie elementary partonic one, 
Eq. flllSp . one should sum up the contributions from quarks and antiquarks weighted 
with the corresponding parton distribution functions (PDFs) in the projectile hadron. 
The hadronic cross-section then reads [OH [96] 

da _ xi da ^ ^ f riA ^ /xi 2\ \ dajgp -» q-f*p) 

3ttM^ xi + X2 Jxi a Va' / d\na 
The PDFs of the projectile enter in a combination which can be written in terms of 
the proton structure function . Notice that with our definitions the fractional quark 
charge Zf is not included in the LC wave function of Eq. (11141) . and that the factor 
in Eq. (11181) accounts for the decay of the photon into the lepton pair. We use the 
standard notation for the kinematical variables Xi and X2 defined in Eq. (llOSp . 

The transverse momentum pt distribution of photon bremsstrahlung in quark- 
nucleon interactions, integrated over the final quark transverse momentum, was derived 
in Ref. P]: 



(119) 



inj L,T' 



where 

S^(x, fi, r2, a) = ^{(Tqq{x, avi) + agg{x, ar2) - (Tqq{x, a(fi - f2))}. (120) 

and ri and ^2 are the quark-photon transverse separations in the two radiation 
amplitudes contributing to the cross-section, Eq. (IllQp . which correspondingly contains 
double-Fourier transformations. The parameter a is the relative fraction of the quark 
momentum carried by the photon, and is the same in both amplitudes, since the 
interaction does not change the sharing of longitudinal momentum. The transverse 
displacement between the initial and final quarks is ari and ar2, respectively. After 
integrating the above equation (I120p over pt, one recovers Eq. (11130 . as one should. 

The hadronic cross-section can then be obtained in the same fashion as given 
in Eq. (IllSp by convolution with the proton structure function. Next we calculate 
the inclusive direct photon spectra within the same framework. For direct photon we 
have M = 0; the transverse momentum distribution of direct photons production from 
hadron-hadron collision reads 

da^ipp^jX) _ 1 da^p^xi ^^ da'^^jq qj) ^^^^^ 



dxpd'^PT Xi + X2 a a d{lna)d'^pT 

We also need to identify the scale Q entering in the proton structure function 
in Eqs. (1118111210 . and relate the variable x of the dipole cross-section entered in 
Eqs. (I113lll20p to measurable variables. From our previous definition, and following 
previous works [HSl EZj, we have that At zero transverse momentum, the 




Figure 37. Left: Inclusive direct plioton spectra obtained ft-om tlie GBW-DGLAP 
dipole models at CDF and CERN energies. Right: The dilepton spectrum with 800- 
GeV beam energy in pp collisions from the E866 |100| fixed target experiment. We show 
the result of the GBW dipole model (dashed hue) and the GBW-DGLAP model (dotted 
line). We also show the result when a constant primordial momentum (A;§) = 0.4GeV^ 
is incorporated within the GBW-DGLAP dipole model (solid line). Experimental data 
are from Refs. [1001 [M]. The figures are taken from [TMIITU5] . 



dominant term in the LC wavefunction Eq. (11141) is the one that contains the modified 
Bessel function Ki{rir). This function decays exponentially at large values of the 
argument, so that the mean distances which numerically contribute are of order l/rj. On 
the other hand, the minimal value of aisxi, and therefore the virtuality which enters 
into the problem at zero transverse momentum is ~ (1 — xi)M^. Thus the hard scale at 
which the projectile parton distribution is probed turns out to be = Pt+ (1 — xi)M^. 
Notice that in the previous studies, [9^ and (1 — xi)M^ [99] were used for the scale 
Q^. Nevertheless, these different choices for Q"^ account for less than a 20% effect at 
small X2 values. 

As example, in Fig. [37] we show the dilepton and inclusive direct photon spectra 
for different experiments. For the dipole cross-section, we use two parametrizations 
introduced in Section 10. It is remarkable that both direct photon production and 
DY dilepton pair production processes can be described within the same colour dipole 
approach without any free parameters. From this study, it is seen that the colour dipole 
formulation coupled to the DGLAP evolution provides a better description of data at 
large transverse momentum compared to the GBW dipole model. 

The colour dipole predictions for the direct photons at the LHC is given in Ref. [103j . 
In the same framework the azimuthal asymmetry of the prompt photons was computed 
in Ref. [104j . for the predictions of other approaches at the LHC see Ref. |105] . 
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Figure 38. Born approximation for elastic scattering in abelian (left) and nonabelian 
(right) theories. 



12. Diffraction 

Diffraction is associated witli tlie optical analogy, which is elastic scattering of light 
caused by absorption. A new feature of diffraction in quantum mechanics is the 
possibility of inelastic diffraction, which is nearly elastic scattering with the excitation 
of one or both colliding hadrons to effective masses which are much smaller that the 
cm. energy of the collision. The main bulk of diffractive events originate from soft 
interactions. Therefore it is still a challenge to describe these processes starting from 
the first principles of QCD. Unavoidably, one faces the problem of confinement which 
is still a challenge for the theory. Nevertheless, the ideas of QCD help to develop quite 
an effective phenomenology for diffractive reactions, i.e., to establish relations between 
different observables. 

12.1. Diffraction in nonabelian theories 

Elastic and inelastic diffraction are large rapidity gap (LRG) processes. Since they 
emerge as a shadow of inelastic interactions, their amplitudes are nearly imaginary. 
This observation is a direct evidence for the underlying theory to be nonabelian. 

Indeed, the elastic amplitude can be mediated only by a neutral exchange in t 
channel, therefore the Born graphs in the abelian and nonabelian cases look like as 
shown in Fig. [38l 

The striking difference between these two amplitudes is in their phases. In the 
abelian case (e.g. in QED) the Born amplitude is real, while in the nonabelian theory 
(QCD) the amplitude is imaginary. 

Data for elastic hadron scattering show that the real part of the elastic amplitude 
is small, and this is a direct evidence for the nonabelian underlying dynamics. This is 
a remarkable observation, since we have so far very few manifestations of nonabelian 
features in the data. 

The Born amplitude depicted in Fig. [38] is independent of energy. Gluon radiation 
gives rise to the energy dependence of the total cross-section through the unitarity 
relation illustrated in Fig. [391 



Applied High Energy QCD 



2Imfgi = gWf% g 



Figure 39. The unitarity relation for the Pomeron amplitude in terms of perturbative 
QCD 

Elastic scattering reaches maximal strength at the unitarity limit of black disc, 
Im/eK&) = l, 

aei = ain = T^R\ (122) 

where R is the radius of interaction. The unitarity relation tells us that the imaginary 
part of the partial amplitude Im/ei(6) cannot rise for ever. After the unitarity bound 
is reached, the total cross-section can rise only due to an energy dependence of the 
interaction radius R{s). Froissart theorem imposes a restriction on this, the interaction 
radius cannot rise with energy faster than R oc ln(s). Then, the total and elastic cross- 
section rise with energy as oc \t?{s) in the Froissart regime of unitarity saturation. 




12.2. Quantum mechanics of diffraction 

Diffractive excitation is a nontrivial consequence of the presence of quantum fluctuations 
in hadrons. In classical mechanics only elastic scattering is possible. An example is 
diffractive scattering of electromagnetic waves. 

One can understand the origin of diffractive excitation in terms of elastic diffraction 
|106[ 1107] . Since a hadron has a composite structure, different hadronic constituents 
interact differently causing a modification of the projectile coherent superposition of 
states. Such a modified wave packet is no longer orthogonal to other hadrons different 
from the incoming one. This makes possible the production of new hadrons, i.e., 
diffractive excitations. 

To simplify the picture, one can switch to the basis of eigenstates of interaction. 
Since a hadron can be excited, it cannot be an eigenstate of interaction, and can be 
expanded over the complete set of eigen states |a) |108l llOQj IllOj : 

\h) = Y.Cl\a), (123) 

a=l 

which satisfy the condition, fei\oi) = /a |«), where fei is the elastic amplitude operator. 

Owing to completeness and orthogonality of each set of states, the coefficient 
in Eq. (11231) satisfy the relations, 

a=l 

= = ^o.p. (124) 

h' 
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The elastic and single diffraction amplitudes can be thus expressed via the eigen 
amplitudes as, 



fh^h _ \ ^ \r'h\2 f 
Jel — I'-^al Jai 



lei 



a=l 



(125) 



a=l 



Using these expressions and the completeness relations, Eqs. fll24p . one can calculate 
the forward single diffraction cross-section without knowledge of the properties of \h'), 
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(126) 



Thus the forward diffractive cross-section is given by the dispersion of the eigenvalues 
distribution. For some specific distributions the dispersion may be zero. For instance if 
all the eigenamplitudes are equal, or one of them is much larger than others. 

According to Eqs. f ll25p -( fT26l) one can calculate the total and diffractive cross- 
sections on the same footing, provided that the eigenstates their weights |C^P 
and the eigenvalues /„ are known. Notice that the eigenamplitudes /„ are the same for 
different hadronic species \ h). This remarkable property of eigen amplitudes is employed 
later on. 

In the Froissart regime all the partial eigen amplitudes reach the unitarity limit, 
Im/a = 1. Then, according to the completeness conditions. 



hh 
el 



h\2 



a=l 



(127) 



a=l 



Diffraction is impossible within a black disc, but only on its periphery, b ^ R. Since in 
the Froissart regime R oc ln(s), 

atot oc Cei oc In^(s), 

cxsd oc ln(s) , (128) 

i.e. <7sd/o'tot oc l/ln(s). The total and single diffractive cross-sections in terms of the 
colour-dipole cross-section read, 

- E \Ca\'^a = I rfVT|vl>,(rT)|V(rT) = {a^rr)) , (129) 

a=l 
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(130) 



where the eigenvalue of the cross-section for a simplest qq dipole (JqqirT) was already 
introduced in Section 10. 
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12.3. DiffracUve DIS 

The contribution of diffractive quasielastic production of vector mesons is a tiny fraction, 
vanishing as l/Q"^, of the total inclusive DIS cross-section. However, the fraction of 
all diffractive events associated with large rapidity gaps in DIS is large, about 10%, 
and is nearly independent of Q^. This turns out to be the result of a contribution 
of rare soft fluctuations in the hard photon. According to Eq. fl99|) a longitudinally 
asymmetric qq pair with a or 1 — a ^ l/Q"^ has a large hadronic size and experience soft 
diffractive interactions like hadrons. Although the admixture of such soft fluctuations in 
the virtual photon is tiny, that may be compensated by a large interaction cross-section. 
This interplay between the fluctuation probability and the cross-section is illustrated for 
inclusive and diffractive DIS in Table. [1] |lllj . 

Table 1. Interplay between the probabilities of hard and soft fluctuations in a highly 
virtual photon and the cross-section of interaction of these fluctuations. 
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a=soft 
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Hard fluctuations of the photon have large weight, but vanishing as l/Q"^ in the 
cross-section, while soft fluctuations have a small, rv?jQ'^ weight, but interact strongly, 
a ~ l/^«g. The latter factor compensates the smallness of the probability in the case of 
DIS, and over-compensates it for diffraction. 

Thus we conclude that inclusive DIS is semi-hard and semi-soft, and the soft 
component is present at any high . On the other hand, diffractive DIS (called 
sometimes " hard diffraction" ) is predominantly a soft process. This is why its fraction in 
the total DIS cross-section is nearly independent. One can test this picture studying 
the dependence of the diffractive DIS [112] . 

Since diffraction is a source of nuclear shadowing [113J, that also should scale in 
X. Indeed, most of experiment have not found any variation with Q'^ of shadowing in 
DIS on nuclei. Only the NMC experiment |114[ 1115] managed to flnd a weak scaling 
violation which agrees with theoretical expectations |116] . 

Notice that in spite of independence of Q^, both diffraction and shadowing are 
higher twist effects. This is easy to check considering photoproduction of heavy flavors. 
In this case the hard scale is imposed by the heavy quark mass, and diffraction becomes 
a hard process with cross-section vanishing as l/mq. Nuclear shadowing also vanishes 
as l/m^ . 
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The true leading twist diffraction and shadowing are associated with gluon radiation 
considered below. 

12.4- Diffractive Drell-Yan reaction 

The dipole description of the Drell-Yan reaction in many respects is similar to DIS, see 
Sections 10 and 11.2. This is not a surprise, since the two processes are related by QCD 
factorization. 

There is an important difference between DIS and DY reaction. In the inclusive 
DIS cross-section one integrates over < a < 1, this is why this cross-section is always 
a mixture of soft and hard contributions (see Table 1). In the case of DY reaction there 
is a new variable, xi, which is the fraction of the proton momentum carried by the 
dilepton. Since a > Xi, one can enhance the soft part of the DY cross-section selecting 
events with xi —>■ 1. This soft part of the DY process is subject to unitarity corrections 
[117] which are more important than in DIS |118j . 

Another distinction between DIS and DY is suppression of the DY diffractive cross- 
section. Namely, the forward cross-section of diffractive radiation qp llqp is zero 
[119] . Indeed, according to Eq. (11261) the forward diffractive cross-section is given by 
the dispersion of the eigenamplitude distribution. However, in both eigen states \q) and 
|g7*) only quark interacts. So the two eigenamplitudes are equal, and the dispersion is 
zero. 

Nevertheless, in the case of hadronic collisions the diffractive DY cross-section does 
not vanish in the forward direction. In this case the two eigen states are \qq) and \qq'^*) 
(for the sake of simplicity we take a pion). The interacting component of these Fock 
states is the qq dipole, however, it gets a different size after the g or g radiate the photon. 
Then the two Fock states interact differently, and this leads to a nonvanishing forward 
diffraction. Notice that the diffractive cross-section is proportional to the dipole size 

12.5. Diffractive Higgs production 

Detection of Higgs particle is the main challenge of the forthcoming experiments of the 
LHC at CERN. The most difficult problem here is to single out a weak signal from high 
backgrounds. One possible process to study is a double diffractive production of Higgs, 
p + p—*p + H + p, with two large rapidity gaps, as illustrated in Fig. HOlfleft). 

Like other diffractive processes, this reaction is strongly suppressed by the small 
survival probability of the gaps. Namely, initial- and final-state inelastic interactions 
of the colliding protons, can easily cause multiparticle production which will fill the 
gaps. The probability of no-interaction is usually called absorptive corrections, which 
are illustrated in Fig. HOT left) by shaded strip. Recent calculations of the cross-section 
of this reaction [122] led to a rather small cross-section, which, nevertheless, may be 
observed due to smallness of the background which is also suppressed by the absorptive 
corrections. 
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Figure 40. Left: Double difFractive Higgs production ^ iJpp. Right: Cross-section 
of exclusive diffractive Higgs production, pp Hpp^ from intrinsic charm (IC), bottom 
(IB) and top (IT) [I2T] . 



Another possible mechanism for Higgs production could be a direct diffractive 
higgsstrahlung similar to diffractive DY. In both cases the radiated particle does not take 
part in the interaction |120j . However, the Higgs coupling to a quark is proportional 
to the quark mass, therefore the cross-section of higgsstrahlung by light hadrons is 
vanishingly small. 

A larger cross-section may emerge due to the admixture of heavy flavors in light 
hadrons. A corresponding mechanism of exclusive Higgs production, pp — ^ Hpp, due 
to direct coalescence of heavy quarks, QQ —>■ H was proposed in Ref. [121] . In this 
case the Higgs is produced not at the mid rapidities, but in the fragmentation region of 
the proton, at large Feynman xp where backgrounds are very small. The cross-section 
of Higgs production was evaluated assuming 1% of intrinsic charm (IC) |123] and that 
heavier flavors scale as l/rriQ [124] . The results are shown in Fig. HOlfright) as a function 
of the Higgs mass for different intrinsic heavy flavors. The cross-section is small, but 
can be detected by dedicated measurements. 



13. Quark and Gluon shadowing 

One may naively expect that the cross-section for scattering a lepton off a nucleus with 
mass number A must be A times bigger than the cross-section for the lepton-proton 
collision. However, several experiments show that the nuclear DIS cross-section at small 
X << 1 is smaller, 

at^ < Aatf. (131) 

This phenomenon is called shadowing. Shadowing has been investigated by various 
experiments in different kinematics ranges. For a review of the experimental and 
theoretical results, see Ref. [125] . 
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Figure 41. Even when nucleons are well separated in the longitudinal direction in 
the infinite momentum frame, gluons fluctuation at small x overlap. 

A particle thrown on a nuclear target has many possibilities of interaction with 
different bound nucleons. However, the total probability of interactions should not 
exceed 1. Therefore, a probability of each interaction must be reduced which can be 
viewed as a result of shadows produced by the preceding collisions. 

Both the Colour Glass Condensate [;74j and shadowing have the same origin: 
longitudinal overlap of gluon clouds originating from different bound nucleons. This 
is illustrated in Fig. HH Bound nucleons in the nucleus do not overlap much, either 
in the rest frame, or in the infinite momentum frame, since both the nucleon size and 
internucleon spacing are subject to Lorentz contraction. However, gluons carrying a 
small fraction x of the proton momentum have a smaller gammafactor and are less 
compressed in the longitudinal direction. Then, the longitudinal propagation of small-x 
partons is large. They overlap and do talk to each other, i.e. they fuse and reduce 
parton density at small x. The cross-section decreases and this is shadowing. Fig. UT] 
shows how gluonic clouds overlap at small x. 

At small x, nuclear scattering is governed by coherence effects which are better 
understood in the target rest frame. As described above, a virtual photon with virtuality 
and energy v splits into a qq pair with a coherence length 

' ~ WTm[, - ^^^^^ 

where M^^ is the effective mass of the fluctuation, and the factor = (1 + M^^/Q^). 
The usual prescription is that M^^ ~ Q"^ since is the only scale available, which leads 
to P = 1/2. Then, the coherence length can be bigger than nuclear radius at low x. 
This means that gg-pair can experience multiple scatterings off different nucleons within 
coherence length. In the infinite momentum frame this corresponds to the overlap of 
parton clouds of different nucleons which leads to diffusion of gluons and consequently 
a reduction of the gluonic density in nuclei. A more careful analysis, however, shows 
that P, even for quarks, depends on the polarization of the photon [98]. The factor P 
for gluons is about one order of magnitude smaller, see Fig. US Therefore, gluons need 
much smaller x in order to overlap in the longitudinal direction. This simple observation 



Applied High Energy QCD 



55 





1 n 




no - 




u.o - 




7 - 


K 
nj 


0.6 - 


S O 




. 


0.5 - 




0.4 - 




0.3 - 




0.2 - 




0.1 - 




0.0 1- 



0.0001 



0.001 



<PT> 



<p°> 



0.01 



0.1 



Figure 42. Bjorken x dependence of P defined in (|132p . corresponding to the 
coherence length for shadowing of transverse and longitudinal photons and gluon 
shadowing, respectively [55|. Solid and dashed curves correspond to = 4 and 
40 GeV^. The bottom curve represent P for gluons. 



leads to a remarkable prediction that the onset of gluon shadowing occurs at smaller x 
compared to quark shadowing. 

The quark and gluon shadowing can be estimated within a simple model at 
high energy. At high energies, dipoles qq are frozen by Lorentz time dilation during 
propagation through the nucleus. Therefore at very small x, it is possible to write entire 
multiple rescattering which occurs during propagator of the qq with fixed transverse size 
r in a eikonal form |126j . 

QAjx) _ 2_ 

where the nuclear thickness function is defined as an integral of the nuclear density 
along the projectile trajectory, T^(&) = J^^dzp{z). Similar calculations can be carried 
out for gluons |127j. 



9a[x} 



- (e^^'^^'^'M'^^^*))) . (134) 



where the gluon-gluon dipole is related to the quark-antiquark dipole cross-section by 
the Casimir factor (Jgg{r) = ^(Jqq{r). Assuming the gluon-gluon fluctuation of the 
projectile have the same distribution function as for qq, one may conclude that the 
effective absorption cross-section providing shadowing is 9/4 times larger than for a qq 
fluctuation of a photon. Such a simple result cannot be true because of the strong gluon- 
gluon interaction which makes their distribution function quite different . Moreover, the 
spin structure of the gluon-gluon distribution function is also different. It turned out 
that in fact gluon shadowing is weaker [llOj . That is because gluons in the proton are 
located within small spots [16], so they have a little chance to overlap in the transverse 
plane, even in heavy nuclei. If the mean value of quark-gluon separation ro, the mean 
number of other dipoles overlapping with this one is 
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Figure 43. Gluon shadowing, for carbon, cooper and lead sX — A GeV^ (solid) 
and Q2 ^ 40 QeV^ (dashed) ^llOj . 



(n,) = ^7rr2(T^)~0.3. (135) 

This indicates that, even at very small x, gluon shadowing must be quite small, see 
Fig. US From experimental data it is very difficult to extract gluon shadowing. For 
the only existing experimental data NMC |114l 1115] , a leading order analysis failed to 
extract the gluon distribution, and the NLO fit turned out to be quite sensitive to gluons 
|128j . Nevertheless, the results indicate a very weak gluon shadowing. 



14. Cronin Effect and nuclear broadening 

Back in 1973, Cronin's group discovered that nuclei may not only suppress reactions, 
but also enhance them [129] . A considerable enhancement was found for production 
of hadrons with large transverse momentum. This effect is measured by the ration i?, 
of the inclusive differential cross-sections for proton scattering on two different targets 
normalized to the respective atomic numbers A and 5, 

BdapA/d^PT 
AdapB / d^PT 

If there were no any nuclear effect, then we had R{pt) = 1; however, for A > B a. 
suppression is observed experimentally at small pt and an enhancement at intermediate 
Pt, and eventually at very high px the ratio seems to approach R{pt) = 1, see Fig. HH 
Experimental data from RHIC |132l I133j for high-p^ hadrons in gold-gold collisions 
raised again the long standing problem of quantitative understanding the Cronin effect. 
In nucleus-nucleus collisions this effect has to be reliably calculated as a baseline for a 
signal of new physics in heavy ion collisions. The only possibility to test models is to 
make comparisons with available data for pA collisions since in pA collisions no hot and 
dense medium is created. 
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Figure 44. Ratio of the charged pion production cross-sections for tungsten and 
beryUium as a function of the transverse momentum of the produced pions [130j . 
Experimental data are from Ref. |131j . 



Soon after the discovery of the Cronin effect, it was understood that the nuclear 
enhancement is a resuh of muhiple interactions in the nucleus |129lll34j . However, in the 
parton model based on QCD factorization, this should be interpreted as a modification 
of PDFs in the nucleus. In the parton model, inclusive particle production for pA 
collisions can be presented in a factorized form, 

^^ = KY.F^/v® F^/A ® i (U - kl) ® D', , (137) 

where da/di{ij kl) is the pQCD parton-parton cross-section and D'^{z,Q''^) are the 
fragmentation functions of a parton k into a hadron h with a fraction z of the parton 
momentum. The i^-factor simulates the NLO contributions. The proton and nucleus 
parton distribution functions were parametrized as 

Fi/p = fi/pixi, Q^) -J- and Fj/a = TA{h) fj/pixj, Q^) -j- . (138) 

where fi/p(A){x, Q"^) are the parton distribution functions of the proton (nucleus). Isospin 
imbalance was taken into account and nuclear shadowing is included by the HIJING 
parametrization [135] . The results of the calculations |136j are depicted in Fig. HS] 

Partons were assumed to have an intrinsic transverse momentum with an average 
squared value {k'^)pA(Ap) and a Gaussian distribution. At a soft scale one does not resolve 
the gluonic structure of a hadron, but only the valence quarks. The mean transverse 
Fermi momentum of these quarks is small (/cq) ~ ^QCD' higher scale, relevant to 
hard reactions, one can resolve the structure of the valence quarks, i.e. the presence of 
gluons and sea quarks. Since those are located at small separations, tq [16], from the 
valence quark, both have more intensive intrinsic Fermi motion, 

{kl) ~ l/rl (139) 
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Figure 45. The photon and vr production cross-sections from the E706 experiment at 
y/s — 31.6 GeV, compared to /cT-corrected NLO calculations |136| . Bottom: the ratio 
(Data- Theory) /Theory for direct photon production. Theory is the NLO calculations 
with primordial parton momentum (kr)- 



This is obviously bigger than the scale associated with nucleon size due to confinement. 
In the parton model it has been shown that even within the next-to-leading order (NLO) 
pQCD correction, experimental data of heavy-quark pair production [137] . direct photon 
production [136] and DY lepton pair production [138] can only be described if an average 
primordial momentum as large as 1 GeV is included (see also Ref. [139] ). For example 
in Fig. HHl the NLO calculations and data for both direct photons and neutral pion 
production are shown [136] . A primordial momentum (/ct) ~ 1-2 GeV seems to provide 
the best description of data. 

A projectile parton propagating through a nucleus experiences multiple interactions 
increasing its transverse momentum. Then the parton participating in a hard collision 
inside the nucleus has an increased transverse momentum compared to Eq. (11391) . which 
corresponds to the interaction with a free proton, 

v^) = {kl) + Ak^b, (140) 

where Ak'^{b, ^/s) is the nuclear broadening. The nuclear broadening is crucial for 
understanding the Cronin effect. Apparently, the strength of the effect depends on 
the relative values of the two terms in Eq. (I140p . In the limit of a weak primordial 
motion the effect should be strongest, while in the case of (k^) » Ak'^ the effect will 
disappear. One may expect Afc^(6, y/s) to be a function of the number of pp collisions, 
i.e. Ak'^{b, i/i) oc app{y/s)TA{b) , where denotes the nucleon-nucleon inelastic cross- 
section. Different parametrizations exist for Ak'^{b, y/s), thought all seem to be rather 
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Figure 46. The probability of multiple interactions via one gluon exchange for the 
quark in the nucleus. The dashed line shows the unitarity cut. 



ad hoc. Here we present a prescription in the framework of colour dipole approach which 
is free from any arbitrary parameters. 

As we aheady mentioned in the previous section, the coherence length Ic is an 
important quantity to understand the effect of multiple parton rescattering [140j . 
Therefore the underlying mechanisms of the Cronin enhancement should also depend 
on the coherence length. In the case of incoherent hard interaction, the incoming 
projectile and outgoing partons experience multiple soft rescattering leading to a high- 
Pt enhancement. At very small x, we are in the large coherence length lc» Ra regime. 
Such a coherent regime is relevant for hadron production at medium large pt at RHIC, 
and it dominates a large range of pt at LHC energies. In addition, in the latter case 
the Cronin effect is substantially reduced by shadowing. 

In the short-coherence length regime Ic << Ra, one can rely on the factorized 
expression Eq. fll37p corrected for broadening Eq. (11401) . The letter can be computed 
within the dipole approach as propagation of a qq pair through the target nucleus. The 
final parton transverse-momentum distribution dNi/(PkiT is written as [1261 1140] : 



^ e-'^{rl+rl){kl) 



: o-q<,(ri-r2,a;)TA(6) 



(141) 

The first bracket in the above equation represents the contribution of the proton intrinsic 
momentum, while the second bracket takes into account the soft parton rescatterings 
on target nucleons. We use the dipole cross-section aqq introduced in Section 10, fitted 
to DIS data. For a gluon when j = g 'm. Eqs. (I137II141I) . we have aqq agg = \crqq. 

Notice that the simple exponential in Eq. (I14ip should not be confused with the 
Glauber eikonal multiple scattering introduced in the previous section. Thus if one needs 
to establish a relation between the expansion of the exponential in the second bracket of 
Eq. (11411) and the multiple quark interaction, it would be incorrect to think that the n-th 
order term of this expansion corresponds to the probability to have n-fold quark multiple 
scattering (we recall the probability cannot be negative!). The appearance of the dipole 
cross-section in Eq. (I14ip . is the result of a product of the amplitude and the time- 
conjugated one, which describe the quarks with different impact parameters. Clearly, 
the object participating in the scattering is not a qq dipole but rather a single coloured 
quark, see Fig. HHl The above prescription describes the fixed target experiments rather 
well, see Fig. HH 
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In the case of a coherence length Ic » Ra, a hard fluctuation in the incident proton 
containing a high-p^- parton propagates through the whole nucleus and may be freed by 
the interaction. Since multiple interactions in the nucleus supply a larger momentum 
transfer than a nucleon target, they are able to resolve harder fluctuations, i.e., the 
average transverse momentum of produced hadrons increases. In this case broadening 
looks like colour filtering rather than Brownian motion. We employ the light-cone dipole 
formulation in the target rest frame which leads to |140] . 

^If"^ (Pt) = fg/p ®a{gA^ g.g^X) ® . (142) 

We assume that high-pr hadrons originate mainly from radiated gluons at such small 
X. The cross-section of gluon radiation reads |110l 11191 1141j 

da{gA 5-15-2^) 



(PpT dyi 



(143) 



where a = p^{gi)/p+{g) is the momentum fraction of the radiated gluon. The function 
Rg incorporates the shadowing effect which originates from the higher Fock components 
\3g), \4:g), etc., missing in the naive eikonalization [1101 1130] . a^g{r,a) is the dipole 
cross-section for a three-gluon colourless system, where f is the transverse separation 
of the final gluons gi and g2- It can be expressed in terms of the usual qq dipole 
cross-sections, 

^3gir) = l{crg,{r) + cr^giar) + or,-,[(l - a)r]} . (144) 

The variable x in cr^(r, a) and Rg is implicit. The light-cone wave function of the gi~g2 
Fock component of the incoming gluon including the nonperturbative interaction of the 
gluons reads |110] . 

.2 



^gg(r,a) = 5- exp 



a{e* ■ e)(e* ■ r) + (1 - a)(e* ■ e)(e* ■ r) 



2rl 

- a(l-a)(e;* ■e2)(e-f)] (145) 

where rg = 0.3 fm is the parameter characterizing the strength of the nonperturbative 
interaction which was fitted to data on diffractive pp scattering. The product of the 
wave functions is averaged in fll43p over the initial gluon polarization, e, and summed 
over the final ones, 61^2- 

In the upper panel Fig. HT] we show the results for RHIC energy ^/s = 200 GeV. 
In the lower panel we show the prediction for the ratio of pion production rates in pA 
and pp collisions obtained using Eqs. f ll42p -( |T^3l) for mid-rapidity at the energy of LHC 
y/s = 5.5 TeV [130] . It is seen that the inclusion of the shadowing effect (sohd line) leads 
to a reduction of the Cronin effect. Note that this curve approaches to unity from below 
at high Pt- We stress that all phenomenological parameters in the above prescription 
are fixed in reactions different from p — A collisions. Therefore, these results may be 
considered free-parameter predictions. 
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Figure 47. Up panel: Ratio of p — Au to pp cross-sections as function of transverse 
momentum of produced pions at the energy of RHIC -y/s = 200. Down panel: 
Predictions for LHC \/s — 5.5 TeV calculated using Eq. (|142p . The dashed and soHd 
curves correspond to calculations without and with gluon shadowing respectively. The 
theoretical curves are taken from ,130J. Data are for d+Au collisions from the PHENIX 
collaboration [142j . 



15. Summary 

During the last half a century, QCD survived through many experimental tests leading 
to a consensus that this is a correct theory of strong interactions. In the asymptotically 
free region, perturbative QCD has been quite successful and many QCD perturbative 
computational tools and techniques have been developed. This is particularly useful in 
order to have a detailed understanding of backgrounds for the search for signals of new 
physics at the LHC. 

Unfortunately, we still have a rather poor understanding of soft nonperturbative 
physics which is never avoidable. Nevertheless, QCD-based phenomenology is well 
developed. Nowadays we are able to calculate many reaction without having to fit 
to the data that we want to explain. On the other hand, the current phenomenology of 
strong interaction looks far more complicated and messy than the first principles (QCD 
Lagrangian) we started with. 

In these lectures we introduced two different approaches to high-energy QCD 
phenomenology: the parton model and the colour dipole formalism. We discussed the 
relevance of both methods as an efficient way to include the non-perturbative features of 
QCD via fitting to some experimental data and predicting others. In the case of parton 
model one fits the universal parton distributions, which then allow one to predict other 
reactions by combining these PDFs with perturbative calculations. Next-to-leading- 
order corrections and higher- twist effects make this program more difficult. In the case 
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of the dipole approach, the universal phenomenological function is the dipole-proton 
cross-section, which is mainly fitted to DIS data from HERA. This description by default 
includes the higher-order and higher-twist corrections. However, this is expected to work 
only at very small Bjorken x and is not useful at large x where valence quarks dominate 
the PDFs. 

The LHC is expected to become a laboratory for gluo-dynamics, which should 
settle many of the controversies in our understanding of small Bjorken x physics. LHC 
data should bring forth important information on the gluonic structures in the proton. 
The currently observed steep rise of the gluon density is expected to be slowed down by 
saturation. This is still debatable, since even in pp at the Tevatron saturation is reached 
only for central collisions. 

The forthcoming LHC data with nuclear beams will reveal the gluonic structure of 
nuclei. They should resolve the controversy about the magnitude of gluon shadowing. 
The saturation scale in nuclei is expected to reach values of a few GeV, leading to strong 
observable effects. 
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